AD-7  72  508 


ASPECTS  OF  MECHANICAL  BEHAVIOR  OF 
ROCK  UNDER  STATIC  AND  CYCLIC  LOADIN'  G 
PART  A:  MECHANICAL  BEHAVIOR  OR  ROCK 
UNDER  STATIC  LOADING 


Jesus  E  Basas ,  et  al 


Wisconsin  University 


y 


N 

Prepared  for: 

Advanced  Research  Projects  Agency 
Bureau  of  Mining 


April  1973 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


ASPECTS  OF  MECHANICAL  BEHAVIOR  OF  ROCK  UNDER 
STATIC  AND  CYCLIC  LOADING 


PART  A:  MECHANICAL  BEHAVIOR  OF  ROCK  UNDER  STATIC  LOADING 


Final  Technical  Report 
April,  1973 


by 


R.W.  Heins  (Co-Principal  Investigator 
with  B.C.  Haim  son) 


Department  of  Metallurgical  &  Mineral  Engineering 

and  the 

Engineering  Experiment  Station 
College  of  Engineering 
The  University  of  Wisconsin 
Madison,  Wisconsin  53706 


ARPA  Order  No.  15  79,  Amendment  3 
Program  Code  No.  2F10.  Contract  No.  H0220041 
Contract  Period:  March  1972  through  April  1973 
Total  Amount  of  Contract:  $50,000 
Sponsored  by  ARPA  and  Monitored  by  U.S.  Bureau  of  Mines 


Discla  imer: 

The  views  and  conclusions  contained  in  this  document  are  those  of  the 
author  and  should  not  be  interpreted  as  necessarily  representing  the 
official  policies,  either  expressed  or  implied,  of  the  Advanced  Research 
Projects  Agency  or  the  U.S.  Government. 


AP 172;  &o  % 


3200.8  (Att  3  to  Enel  1) 

_  Mar  7»  86 

^fnnly  (‘1  il  r.un.n  _ _ — — — — — — — — — — — — — 

,  DOCUMENT  CONI  KOL  DAI  A ■ H  D 

,  , frrn.lt,  rl..  II.  fan  ■(  IW  .  ■  I  . .  ,.o. ...toll.. . I  t  ml.  .1  ->■<■»■  »■«  .fp»M  I.  rlnf .!»<«) 

■■ ;vjt.,T7.«.”  u.  ,..  .  o.,, ..  CL.5..r.c*no.. 

Engineering  Experiment  Station  Unclassified _ 

University  of  Wisconsin  ,t>'  *"0,,r> 

Madison,  Wisconsin  53706 _ — - - 


j  iirpon  r  t  1 3  i  r 


Aspects  of  Mechanical  Benavior  of  Rock  under  Static  and  Cyclic  Loading 
Part  A  -  Mechanical  Behavior  of  Rock  under  Static  Loading 


4  UtSCKlRTIVl  NO  I  ».  »  (  T>  k>1  mnd  inch'**  dBU’B) 

Annual  Technical  Report  _ _ _ — _ - 


*.  At  i  iiONitl  filfil  tnliiml,  Ittl  n»mi)  , 

Jesus  E.  Basas,  Robert  W.  Heins 


t|»*OM  T  O  *  T  t 

April  19  73 


•  «.  CONTRACT  OR  GRANT  NO 

H022004 1 

6.  p«ojrc r no  ARPA  Order  No.  1579 
Amendment  3 
‘  Program  Code  No.  2F10 


7B.  tOTAL  NO  OF  PACri  NO  OF  RftF* 

99  4 


|  Bm.  ORIGINATOR**  REPORT  NUMRCNlM 

144-C745 


lb.  OThlN  H4.RORT  NO|V>  (Any  0lh»t  WllfcAN  I/Ml  MT  bB  BBBl$nmd 
hi*  tepofl) 


IO  LHtTRlRUTIQN  ITATl^TNT 


Distribution  of  this  document  is  unlimited. 


II  UiCPtlMLNURV  NOTLI 


Ij  AOSTRACT 


1W  IPONIUHINC  MILITARY  ACTIVITY 


Advanced  Research  Projects  Agency 
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A  finite  element  program  based  on  both  2-  agd  3-dimensional  element 
models  is  developed  to  simulate  the  Brazilian  test.  Rock  isotropy,  non- 
homogenetty,  and  anisotropy  are  taken  into  account  in  the  analysis.  The  pro¬ 
gram  has  two  failure  criteria  options:  the  old  'criterion  described  in  the  semi¬ 
annual  report  and  a  new  criterion  in  which  the  elastic  modulus  across  tension 

cracks  is  reduced  to  zero.  Tests  runs  have  shown  that  the  new  failure  criterion 

< 

yields  far  more  realistic  load-displacement  and  load-strain  curves. 

A  new  equation  solving  procedure,  the  sparse  matrix  technique,  is 
used  in  the  program  making  possible  the  solution  of  certain  problems  which 
could  not  be  handled  by  the  banded  matrix  technique  because  of  computer 
storage  limitations.  To  demonstrate  the  capabilities  of  the  program,  severer 
test  problems  were  run.  The  results  are  shown  in  Chapter  2  of  this  report. 
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PREFACE 


This  report  covers  the  second  year  accomplishments  in  the  research  pro¬ 
gram  entitled,  "Mechanical  Behavior  of  Rock  Under  Static  Loading,"  R.W. 

Heins,  Co-Principal  Investigator.  The  program  is  Part  A  of  the  project  entitled, 
"Aspects  of  Mechanical  Behavior  of  Rock  Under  Static  and  Cyclic  Loading" 
(Contract  No.  H022004 1).  Part  B  of  the  project  is  published  in  a  separate  volume. 
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SUMMARY 


ASPECTS  Or  MECHANICAL  BEHAVIOUR  OF  ROCK  UNDER  STATIC  LOADING 

PART  A 

Summary  of  Year's  Work 

Brazilian  tests  were  carried  out  on  dacite,  Valder's  limestone,  and 
St.  C’oud  gray  granodiorite  to  determine  size  effect  on  tensile  strength. 
Curves  showing  variation  of  tensile  strength  with  specimen  dimensions  and 
other  findings  constitute  Chapter  1  of  the  semi-annual  technical  progress 
report^. 

In  connection  with  the  theoretical  phase  of  the  investigation,  a  2- 
dimensional  computer  program  which  simulates  the  Brazilian  test  is  developed. 
The  program  is  based  upon  isoparametric  finite  elements  which  may  be  either 
isotropic  or  non-homogeneous  and  a  failure  criterion  in  which  a  fixed  percent¬ 
age  of  the  stiffnesses  of  the  failed  elements  is  subtracted  from  the  total  stiff¬ 
ness  during  each  loading  cycle.  The  complete  program  listing  and  the  results 
of  several  computer  runs  are  part  of  the  semi-annual  report^1  \ 

A  second  program  is  presented  in  this  report.  The  program  employs 
2-  and  3-dimensional  isoparametric  elements  which  exhibit  anisotropy  as 
well  as  isotropy  and  non-homogeneity.  Two  failure  criteria  are  used:  the 
old  criterion  described  in  the  preceding  paragraph  and  a  new  criterion  in  which 
the  elastic  modulus  across  tension  cracks  is  reduced  to  zero.  Based  on  few 
test  runs  made,  the  new  failure  criterion  yields  considerably  superior  and 
more  realistic  load-displacement  and  load- strain  curves. 

To  demonstrate  the  capabilities  of  the  program,  several  problems  were 
run  and  the  results  are  shown  in  Chapter  2.  The  use  of  a  new  equation  solver, 
the  sparse  matrix  technique,  permitted  the  solution  of  certain  problems  which 
could  not  be  handled  by  the  banded  matrix  technique  because  of  the  bandwidth 


limitation. 
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CHAPTER  1 


PROGRAMMING  ROUTINES 


1 . 1  Introduction 

A  finite  element  program  for  the  prediction  of  rock  fracture  charac¬ 
teristics  in  Brazilian  tests  is  presented  herein.  The  program  has  all  the 
capabilities  of  the  previous  2-dimensional  program  (1)  plus  the  following: 

1.  Availability  of  both  2-  and  3-dimensional  elements 
with  isotropic,  anisotropic,  or  non-homogeneous 
material  property. 

2.  A  failure  criterion  in  which  the  elastic  modulus 
across  tension  cracks  is  reduced  to  zero. 

3.  An  equation  solver  in  which  only  the  non-zero  elements 
of  the  stiffness  matrix  are  involved  in  the  elimination 
and  back-substitution  processes. 

4.  A  nodal  reordering  scheme. 

The  theoretical  considerations  upon  which  the  program  is  based 
have  already  been  discussed  in  the  semi-annual  report  (1)  and  are  not  re¬ 
peated  here.  The  program's  major  routines  are  described  in  the  following 
sections. 

1 . 2  Nodal  Reordering 

Among  the  different  mathematical  undertakings  encountered  in  a  finite 
element  program,  the  solution  of  the  governing  linear  equilibrium  equations 
determines  to  a  large  extent  the  effectiveness  of  the  program.  When  the 
first  finite  element  programs  were  developed  only  small-size  systems  could 
be  treated  because  of  limited  memory  space  for  the  storage  of  the  global 
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stiffness  matrix.  With  the  introduction  of  external  storage  disks,  it  be¬ 
came  possible  to  store  only  a  small  part  of  the  stiffness  matrix  in  com¬ 
puter  memory  at  a  time,  the  balance  being  in  the  disk,  permitting  the  solu¬ 
tion  of  considerably  larger  systems.  The  most  widely  used  disk-aided 
equation  solving  procedure  is  the  banded  matrix  technique  developed  in 
1965  by  \:.L.  Wilson  of  the  University  of  California.  Wilson's  method 
takes  advantage  of  the  banded  property  of  the  stiffness  matrix  resulting  in 
savings  in  execution  time  and  storage  requirement.  Its  only  limitation  is 
in  the  width  of  the  band.  Unfortunately,  the  maximum  bandwidth  that  most 
digital  computers  allow  is  not  large  enough  in  many  practical  situations. 

In  the  context  of  the  present  research,  such  situations  arise  when  the  whole 
circular  face  of  the  Brazilian  test  specimen  is  considered  in  the  analysis 
and/o  when  3-dimensional  problems  are  treated.  In  both  of  these  cases, 
because  of  the  closed  geometry  of  the  finite  element  mesh  and  the  high 
nodal  connectivity  of  the  3-dimensional  elements,  a  low-numbered  node  be¬ 
comes  connected  to  a  high-numbered  node  resulting  in  prohibitively  large 
bandwidths.  Nodal  reordering  schemes  have  been  developed  to  minimize 
the  bandwidth.  The  effectiveness  of  such  schemes,  however,  is  very  limited. 

An  equation  solver  which  stores  and  processes  only  the  non-zero 
elements  of  the  stiffness  matrix  is  employed  in  the  present  program.  To 
minimize  the  number  of  non-zero  elements  created  during  the  elimination 
process,  a  reordering  of  the  nodes  is  undertaken.  The  reordering  concept 
used  is  based  upon  one  described  by  Jensen  and  Parks  (2).  The  complete  re¬ 
ordering  flow  chart  is  shown  in  Figure  1.3.  To  start  with,  the  nodes  are  as¬ 
signed  numbers  in  the  manner  shown  in  Figures  1.1  and  1.2.  This  allows  for 
the  automatic  generation,  within  the  program,  of  the  nodes  and  their  con¬ 
nections .  The  reordering  of  the  nodes  then  starts.  At  each  stage  of  the  re¬ 
ordering  process,  the  branch  table  of  the  last  node  (JL0C)  entered  in  the 


Figure  1.2  Numbering  Scheme  for  Nodes  in  Three-dimensional  Mesh 


Figure  1.3.  Flow  Chart  of  Nodal  Reordering  Routine. 
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fp ordering  array*  is  scanned  for  the  node  with  the  fewest  number  of 
connections.  If  such  node  has  equal  or  fewer  connections  than  JL0C, 
then  it  is  placed  into  the  next  available  space  of  the  reordering  array. 
Otherwise,  the  node  with  the  fewest  number  of  connections  from  among 
those  which  remain  to  be  renumbered  is  placed  into  the  next  available 
space  of  the  reordering  array.  The  process  is  repeated  until  all  nodes 
have  been  renumbered.  Whenever  more  than  one  node  satisfy  the  afore¬ 
mentioned  conditions,  the  node  with  the  lowest  number  is  selected.  An 
example  of  the  effect  of  reordering  is  illustrated  in  Figures  1.4  and  1.5. 

In  Figure  1.4c  no  reordering  is  undertaken  and  18  non-zero  elements 
(represented  by  the  letter  c)  are  created  during  the  elimination  process. 

In  the  reordered  version.  Figure  1.5c,  only  11  such  elements  are  formed. 

The  effect  of  reordering  becomes  even  more  significant  in  large  systems 
v/here  a  high  percentage  of  the  non-zero  elements  are  created  during 
elimination. 

The  removal  of  JL0C  from  the  branch  tables  of  the  nodes  connected 
to  it  (Step  8  of  flow  chart)  is  equivalent  to  the  elimination  of  non-zero  ele¬ 
ments  located  directly  below  the  diagonal  element  in  the  row  of  JL0C.  Nodes 
added  to  the  same  branch  tables  (Step  10)  on  the  other  hand,  represent  the 
non-zero  elements  created  as  a  result  of  the  activity  in  Step  8.  Steps  8 
and  10  make  up  the  pseudo-elimination  process  which  is  an  essential  part 
of  the  reordering  routine.  The  final  contents  of  the  branch  tables  represent 
the  off-diagonal  non-zero  elements  of  the  upper  triangular  form  of  the  stiff¬ 
ness  matrix.  The  location  of  these  elements  has  to  be  defined  before  the 
actual  elimination  can  be  carried  out. 

Two  arrays  are  needed  to  locate  a  non-zero  element.  These  are  NR0W, 
generated  for  convenience  during  the  actual  elimination  process  (Figure  1.6), 


*  NRL0  in  which  NRC0(i)  gives  the  original  number  of  the  node  assigned 
the  new  number  i. 


(a)  Finite  Element  Mesh 
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Figure  1.4.  Unreordered  Stiffness  Matrix 
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Figure  1.5.  Reordered  Stiffness  Matrix 
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and  L0C,  generated  in  the  reordering  routine.  NR0W(j)  gives  the  number 
assigned  to  the  first  off-diagonal  non-zero  element  in  row  j.  L0C(i)  gives 
the  column  location  of  the  ith  non-zero  element.  Thus  the  kth  non-zero 
element  would  be  located  in  column  L0C(k)  and  in  row  m  where  NR0W(m)  < 
k  <  NR0W (m+1) .  Sequential  numbers,  starting  with  one,  are  used  to  identify 
the  non-zero  elements.  The  numbering  proceeds  in  the  row -wise  direction 
with  the  off-diagonal  elements  in  a  row  nunbered  ahead  of  the  diagonal  ele¬ 
ment.  The  following  tables,  used  to  locate  the  non-zero  elements  in 
figure  1.5c,  serve  to  clarify  the  idea  presented  in  this  paragraph.  The  ar¬ 
rows  point  to  the  diagonal  elements. 

*  *  *  * 
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Table  1.1  Locator  Arrays  for  Non-zero  Elements  of  Figure  1.5c 
1.3  Gaussian  Elimination 

In  the  discussion  in  the  preceding  section,  it  should  be  noted  that  a 
'non-zero  element"  is  actually  a  submatrix  whose  order  is  equal  to  the 
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number  of  degrees  of  freedom  specified  at  a  node  and  a  "row"  are  actually 
rows  occupied  by  the  submatrices.  Hereafter,  these  expressions  shall 
convey  the  same  meanings.  Thus,  a  "non-zero  element"  is  either  a  2x2 
submatrix  (2-dimensional  cases)  or  a  3x3  submatrix  (3-dimensional  cases). 
Similarly,  the  "nth  row"  of  the  stiffness  matrix  would  be  taken  to  mean  the 
nth  2  rows  or  the  nth  3  rows  of  the  matrix. 

The  flow  chart  of  the  Gaussian  elimination  routine  is  shown  in 
Figure  1.6.  In  Step  3,  only  the  non-zero  elements  on  and  to  the  right  of 
the  diagonal  need  be  stored  in  £  s]  because  of  symmetry.  The  elements  are 
stored  in  the  same  order  that  they  would  appear  in  the  stiffness  matrix  with 
enough  spaces  reserved  for  the  non-zero  elements  which  will  be  created 
later  on  during  the  elimination  process.  Using  Figure  1.5c  again  as  an  ex¬ 
ample,  [s]  would  appear  as  follows  during  the  storage  of  the  different  rows 
of  the  stiffness  matrix: 


Row  1: 


11  S 12 

S 13  S14  S 15 

S 16  S17  S 18 

21  S22 

S23  S24  S25 

S26  S27  S28 

Row  6: 


Row  9: 


S11 

S 1 2 

0  0 

S15 

S16  S17 

S18  S19 

Sl,  10 

S21 

S22 

0  0 

S25 

S26  S27 

S28  S29 

S2,  10 

S11 

S12 

0  0 

0  0 

S17  S18 

0  o“ 

S21 

S22 

0  0 

0  0 

S27  S28 

0  0 

13 


Row  16: 


in  which  the  spaces  marked  s  are  occupied  by  the  x  elements  while 
those  marked  0  (zero)  are  reserved  for  the  c  elements. 

After  the  formation  of  any  row,  say  i,  the  process  of  eliminating  all 
non-zero  elements  in  the  row  located  to  the  left  of  the  diagonal  immediately 
follows.  Had  the  entire  upper  triangular  half  of  the  stiffness  matrix  been 
considered  in  the  elimination  procedure,  these  elements  would  Le  simply  the 
"mirror  image"  of  the  elements  lying  on  the  ith  column  of  the  preceding  rows. 
Since  this  is  not  the  case,  a  method  for  locating  the  aforementioned  elements 
has  to  be  devised.  This  is  accomplished  by  introducing  the  spotter  array 
NSPT.  The  spotter  array  keeps  track,  for  each  row,  of  the  element  in  that 
row  whose  "mirror  image"  is  next  in  line  for  elimination.  Because  of  the 
row-wise  direction  of  the  elimination  procedure,  the  order  in  which  the  ele¬ 
ments  of  any  particular  row  are  used  in  the  elimination  process  proceeds  from 
left  to  right;  that  is,  the  "mirror  image"  of  the  left-most  element  is  eliminated 
first  and  that  of  the  right-most  element  is  eliminated  last.  Thus,  at  the  start 
of  the  elimination  process  NSPT(j)  =  NRCfW(j)  for  any  row  j  and  after  each  elim¬ 
ination  the  NSPT  value  of  the  affected  rows  is  incremented  by  one.  NR0W(j) 
is  defined  on  page  11.  The  procedure  discussed  in  this  paragraph  corresponds 
to  Steps  4  through  6  of  the  flow  chart  shown  in  Figure  1.6. 

It  should  be  noted  that  after  elim  nation  of  a  row,  the  elements  in  that 
row  are  stored  in  the  next  available  space  in  tape  or  drum  #10.  Normalization 
of  an  element  (pre -multiplication  of  the  element  by  the  inverse  of  the  diagonal 
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(4) 


(5) 


(6) 

(7) 


(8) 


(9) 

(10) 


Figure  1.6  Continued. 
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element)  takes  place  only  when  the  element  is  no  longer  needed  in  the 
elimination  process  (Step  9).  The  actual  row  operations  carried  out  during 
the  elimination  process  are  represented  by  Step  8. 


The  graph  of  F(S),  with  the  maxima,  minima,  and  inflection  points  indicated, 
is  shown  in  Figure  1.7. 

Applying  the  Newton-Raphson  method,  the  roots  of  equation  (1.1)  can 
be  evaluated  by  means  of  the  iteration  formula 


F(sn) 

&n+l  Sn  F'(s  )  U.2) 

n 

in  which 

F'(S)  =  df"  =  3s2  “  2aiS  +  a2 

n  =  0,  1,2 . 

S^  =  approximation  of  root  value  at  nth  iteration 


Setting  Sq  -  —  (that  is,  starting  the  iteration  at  the  inflection  point  I), 

equation  (1.2)  will  yield  the  value  of  the  middle  root  R  With  R  known, 

2  2 


equation  (1. 1)  may  be  rewritten  as 
F(S)  =  (S-R2)Q(S)  =  0 


(1.3) 


in  which  Q(S)  is  a  quadratic  polynomial  whose  roots  are  the  two  remaining 
roots  of  equation  (1.1).  Q(S)  may  be  conveniently  written  as 


Q  (S)  =  S*  +  bjS  +  b2 


(1.4) 


in  which  (4) 


bl  =  R2  -al 


b2  R2bl+a2 


Thus  the  two  remaining  roots  are 


Ri  ‘ 


-bjVbj 


R3  ’ 


-bj+Vbj  -4b2 


(1.5) 


The  direction  cosines  1  m  n  of  the  principal  stresses  are  ob- 

1  1  #  1 

tained  from  the  equations 


(Ri  -  V'i  -  =  0 


TrS1i  +  (R|-  Vmi-  ezni  * 


0 


0 


(1.6) 


"rz*i  tQz 


m , 


(R.  - 
1 


°z)ni  = 


taking  into  account  the  relationship 

2  2  2 

if  +  mf  +  n  =  1 

i  i  i 


(1.7) 


MM 
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CHAPTER  2 

NUMERICAL  EXAMPLES 


2 . 1  Introduction 

To  demonstrate  the  capabilities  of  the  program,  several  cylinders 
exhibiting  all  three  rock  properties  considered  in  the  program  are  analyzed. 
Unfortunately,  it  is  not  possible  to  compare  the  results,  particularly  with 
respect  to  the  three-dimensional,  anisotropic,  and  non-homogeneous  cases, 
with  any  existing  alternative  solution,  so  relatively  coarse  meshes  are 
used  to  save  on  computer  expense  which  quite  prohibitive  in  most  cases. 
Isotropic  problems  are  run  mainly  to  compare  the  load-displacement  and 
load -strain  curves  ootained  by  the  two  failure  criteria  used  in  the  program. 

In  the  new  failure  criterion,  the  stiffness  matrix  of  the  elements  that  failed 
is  revised  by  reducing  to  zero  the  elastic  modulus  across  tension  cracks 
when  failure  is  in  tension  or  by  removing  a  fixed  percentage  (specified  by 
program  user)  of  said  matrix  when  the  failure  is  in  compression  or  shear. 

In  the  old  failure  criterion  (used  in  semi-annual  report) ,  the  stiffness  matrix 
of  the  failed  elements  is  reduced  by  a  fixed  percentage  regardless  of  whether 
tensile,  compressive  or  shearing  failure  occurs.  One  shortcoming  of  the  old 
criterion  is  the  unrealistic  shape  of  the  resulting  load-strain  curves  and,  to 
a  certain  extent,  load-displacement  curves. 

In  anisotropic  discs  and  cylinders  the  allowable  stresses  at  a  point 
vary  with  direction  making  the  determination  of  the  load  factor  (critical  ratio 
of  actual  stress  to  allowable  stress)  of  an  element  quite  involved,  one  re¬ 
quiring  time-consuming  trial-and-error  procedure.  Because  of  this,  a  simplify¬ 
ing  assumption  is  made  in  the  program  in  which  the  critical  direction  is  con¬ 
sidered  to  coincide  with  one  of  the  principal  stresses  as  in  the  case  of  iso¬ 
tropic  discs  and  cylinders.  This  assumption  makes  it  necessary  only  to 


A 


determine  the  allowable  stress  in  the  direction  of  the  critical  principal 
stress  to  obtain  the  element  load  facfor.  In  the  program  the  allowable 
stresses  at  a  point  are  assumed  to  vary  as  an  ellipsoidal  surface  with  the 
principal  axes  coinciding  with  the  directions  of  anistropy. 

<  onsiderable  effort  was  expended  to  check  each  individual  subroutine 
of  the  program  by  long  hand  calculation.  As  part  of  the  debugging  process, 
it  was  also  decided  to  rerun  thi  isotropic  problem  presented  as  Problem  2 
in  the  semi-annual  report.  The  results  were  in  close  agreement  with  those 
obtained  by  the  previous  program. 

All  the  cylinders  treated  herein  possess  the  following  dimensions: 

Diameter  =  3  inches 
Length  =  1  inch 

"Vertical  displacement"  refers  to  the  average  displacement  of  the  loaded 
nodes  in  the  direction  of  the  load  while  "horizontal  strain"  refers  to  the 
strain  at  the  central  point  of  the  disc  or  cylinder  measured  normal  to  the 
loaded  plane. 
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2 . 2  Two-dimensional  Anisotropic  Problem 


Elastic  Moduli: 

E^  =  3,000,000  psi 

E^  =  5,700,000  psi 

Allowable  Compressive  Stresses: 

C  =  15,000  psi 
C2  =  27,000  psi 

Allow.  Tension  =  .05  of  Allow.  Compression 
Allow.  Shear  =  .  10  of  Allow.  Compression 

Poisson's  Ratio  =  .25 
Shear  Modulus  =  1,600,000.  psi 

Finite  Element  Mesh  (Whole  Disc): 

5  radial  divisions 

16  circumferential  divisions 

Directions  of  anisotropy  are  indicated  in  Figure  2.1. 

2  • 3  Two-dimension  Isotropic  Problem  (Old  Failure  Criterion) 

Elastic  Modulus  =  5,700,000.  psi 
Allow.  Compressive  Stress  =  27,000.  psi 
Allow.  Tens!  n  =  .05  of  Allow.  Compression 
Allow.  Shear  =  .10  of  Allow.  Compression 
Poisson's  Ratio  =  .25 
Finite  Element  Mesh  (Disc  Quadrant): 

6  radial  divisions 

9  circumferential  divisions 

Stiffnefs  matrices  of  failed  elements  are  reduced  90% 


INPUT  DATA  LISTING  FOR  PROBLEM  2 


Critical  Load,  lbs. 


Vertical  Displacement,  10_^  inch 
Horizontal  Strain,  10“ 5  inch/inch 


Figure  2.4.  Load-Displacement-Strain  Curves  (Problem  2.3) 
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2 . 4  Two-dimensional  Isotropic  Problem  (New  Failure  Criterion) 

Data  are  the  same  as  Problem  2.3.  Stiffness  matrices  of  elements  failing 

in  compression  or  shear  are  reduced  90%.  Elastic  moduli  across  tension  cracks 
are  reduced  to  zero. 

2.5  Three-dimensional  Isotropic  Problem  (Old  Failure  Criterion) 

Elastic  Modulus  =  5,700,000.  psi 
Allow.  Compressive  Stress  =  27,000.  psi 
Allow.  Tension  =  .05  of  Allow.  Compression 
Allow.  Shear  =  9.00  of  Allow.  Compression 
Poisson's  Ratio  =  .25 
Finite  Element  Mesh  (Cylinder  Quadrant): 

5  radial  divisions 

6  circumferential  divisions 
2  axial  (Z-axis)  divisions 

Because  of  symmetry  only  half  the  length  of  cylinder  is  considered 
in  analysis.  Stiffness  matrices  of  failed  elements  are  reduced  90%. 


2.6  Three-dimensional  Isotropic  Problem  (New  Failure  Criterion) 

Data  are  the  same  as  Problem  2.5.  Stiffness  matrices  of  elements 
failing  in  compression  or  shear  are  reduced  90%.  Elastic  moduli  across  ten¬ 
sion  cracks  are  reduced  to  zero. 


A 


i 

i 
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P 


♦Failed  In  shear.  Rest  failed  in  tension. 
Figure  2.5.  Progression  of  Failure  In  Problem 


Note:  All  elements  failed  in  tension. 


Figure  2.7.  Progression  of  Failure  in  Problem  2.5. 


Critical  Load,  lb/inch 


+ - ♦ - • - 1 - 1 - I - ( - 1 - » 

*234  5  6789 

Horizontal  Strain,  10“5  inch/inch 


Figure  2.8.  Load-Displacement-Strain  Curves  (Problem  2.5). 


INPUT  DATA  LISTING  FOR  PROBLEM  2 


2.7  Three-dimensional  Anisotropic  Problem 


Elastic  Moduli: 

E^  =  3,000,000.  psi 

E2  =  5 , 700,000.  psi 

E3  =  5 , 700, 000.  psi 

Allow.  Compressive  Stresses: 

C1  =  15,000.  psi 

C2  =  27, 000.  psi 

=  27, 000.  psi 

Allow.  Tension  =  .05  of  Allow.  Compression 

Allow.  Shear  =  .10  of  Allow.  Compression 

Poisson's  Ratios: 

/i-L  =  .25 

/i2  =  .23 

/i3  =  .23 

Shear  Moduli: 

G ^  =  1, 600, 000  psi 

G2  =  2,800,000  psi 

G^  =  2 , 800, 000  psi 

Finite  Element  Mesh  (Whole  Cylinder): 

3  radial  divisions 
8  circumferential  divisions 
2  axial  (Z-axis)  divisions 

Directions  of  anisotropy  are  indicated  in  Figure  2.11. 


INPUT  DATA  LISTING  FOR  PROBLEM 
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Front  element  failed. 

Both  front  and  back  elements 
failed . 


3  =  °° 

Directions  of  Anisotropy 


Figure  2.11.  Progression  of  Failure  in  Problem  2.7. 
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0  12  3  4 

Horizontal  Strain,  10~5  lnch/lnch 


Figure  2.12.  Load-Displacement-Strain  Curves  (Problem  2.7) 


2.8  Three-dimensional  Non -homogeneous  Problem 

Elastic  Modulus  =  5 ,500, 000  to  6. 500, 000  psi 

Allow.  Compressive  Stress  =  22,000  to  32,000  psi 

Allow.  Tension  =  .05  of  Allow.  Compression 

Allow.  Shear  =  .  10  of  Allow.  Compression 

Poisson' s  Ratio  =  .23  to  .25 

Finite  Element  Mesh  (Whole  Cylinder): 

3  radial  divisions 
8  circumferential  divisions 
2  axial  divisions 
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I 


Figure  2.13,  Progression  of  Failure  in  Problem  2.8. 
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2 . 9  Discussion  of  Results  and  Conclusion 

It  is  difficult  to  verify  explicitly  the  theoretical  correctness  of  the 
present  solution  due  to  lack  of  an  alternative  solution  to  compare  the  re¬ 
sults  with.  To  compensate  for  this  situation,  extensive  long-hand  checks 
were  undertaken  on  the  individual  routines  that  make  up  the  program.  The 
results,  on  the  whole,  are  fairly  reasonable  and  seem  to  indicate  the  cor¬ 
rectness  of  the  whole  program  and  the  theory  upon  which  it  is  based. 

The  unrealistic  failure  shown  in  Fig.  2.1  could  be  attributed  to  the 
following: 

1)  Coarseness  of  the  finite  element  mesh. 

2)  The  simplifying  assumption  used  in  anisotropic  cases  (second 
paragraph  of  section  2.1)  that  the  critical  direction  is  neces¬ 
sarily  along  one  of  the  principal  stresses. 

Item  1)  above  could  also  be  the  cause  of  the  unrealistic  failure  patterns 
shown  in  Figs.  2.3  and  2.5,  It  should  be  noted  that  although  the  meshes 
shown  in  these  figures  appear  to  be  relatively  fine,  they  are  actually  approxi¬ 
mately  six  times  coarser  than  the  mesh  of  Problems  two  and  three  of  the  semi¬ 
annual  report  (see  page  49  of  same  report).  Unrealistic  failure  pattern  were 
also  observed  during  the  debugging  of  the  previous  two-dimensional  program 
when  coarse  meshes  were  used. 

Problems  2.3  through  2.6  were  run  to  determine  the  relative  merits  of 
the  two  failure  criteria  used  in  the  program.  There  is  no  significant  dif¬ 
ference  between  the  two  as  far  as  the  failure  patterns  are  concerned.  The 
displacement  and  strain  curves,  however,  differ  considerably  particularly 
in  the  coarser  three-dimensional  mesh  (see  Figures  2.8  and  2.10).  The  new 
failure  criterion  seems  to  yield  consistently  definable  curves  while  the  old 
criterion  does  not.  Thus  although  the  failure  patterns  are  unrealistic,  the 
plots  show  a  much  better  relationship  between  load  and  strain  and  load  and 
displacement. 


Table  2. 1  shows  the  results  of  the  nodal  reordering  carried  out  in 
connection  with  Problem  2.7.  It  took  no  more  than  30  seconds  to  complete 
the  whole  reordering  process. 

Table  2.2  gives  some  information  concerning  the  sizes  of  the  test 
problems  described  in  this  chapter.  The  maximum  number  of  nodes  and 
non-zero  elements  of  the  stiffness  matrix  that  the  program  can  handle  are 
indicated  in  Table  2.3.  These  limits  may  still  be  increased  depending  cn 
the  computer  used. 

All  mathematical  operations  in  the  program  are  carried  out  in  single 
precision.  It  may  be  necessary  to  resort  to  double  precision  to  reduce 
rounding-off  errors. 


APPENDIX  A 
PROGRAM  LISTING 

The  program  listed  in  the  following  pages  is  written  in  Fortran  V 
language  for  use  specifically  in  the  UNIVAC  1108  Computer  located  at  the 
University  of  Wisconsin  -  Madison  campus.  It  requires  six  external  storage 
units  (either  tape  or  drum)  which  are  assigned  the  numbers  10,  11,  12,  13, 
14,  and  15.  The  contents  of  the  files  written  in  units  10  and  12  continually 
undergo  changes  during  a  program  run  making  it  necessary  for  these  particu¬ 
lar  units  to  posses  random  access  capabilities. 

In  addition  to  the  subroutines  listed,  two  other  subroutines, 

RANUN(R)  and  RANUNS(N),  are  called  in  the  program  during  execution  of  non- 
homogeneous  problems.  RANUN(R)  generates  pseudo-random  numbers  as¬ 
suming  a  statistically  uniform  distribution  while  RANUNS(N)  sets  to  N  the 
starting  point  of  the  random  number  generator;  N  is  supplied  by  the  user. 

Both  subroutines  were  developed  by  the  Madison  Academic  Computing  Center 
(MACC)  and  are  part  of  the  MACC  Random  Number  Routine  package. 

Except  for  the  replacement  of  RAN  UN  (R)  and  RANUNS(N),  there  should 
be  no  difficulty  in  making  this  program  adaptable  to  other  computer  systems. 
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MAIN  PROGRAM 


COMMON/rUNKY/NSRCjSOO.NRE  0(500)  .N0ER1500) 
COHNON/RUTH/kOC  (  25000 )  . HTnh.M^NIUM  INUM 
COMMON/ JNrER /NR, NS, NZ ,NWOP,NDIM 
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133. 

134. 
13ft. 

114. 
137. 
Lift. 
13f. 
no, 
m. 
112. 

ns, 
n<t, 
n». 
n*. 

n7, 

ne_a 

nt. 
1*0. 
ifti, 
1*2. 
113. 

_!»«• 

IS*. 

1*6. 

187. 

LI*. 

1*4. 

140, 

141, 

142, 

143, 

144, 
14ft. 
144. 
144, 

148, 

144, 


PRINT  313 
PRINT  31N.01A 
PRINT  315, H 
PRINT  3 1 4  ,NELEM 
PRINT  3?0  *  NTNN 

PRINT  SlA.lRCm  ,Ifl,NRl  »  _ _ 

PRINT  Jl7,lTCU).nl,Nf  '  I 
DO  15  1  «2  ,NS 1 

IS  TCItmei  1  >•. 0174532428 
1 f  <  NO  l M  .  EQ  i  2)60  TO  14 
PRINT  318,(20(1) .!*1,NZ1  ) 

14  CALL  NDTNDIiACilNGA iLNfiO )  _ 

CALL  RCOROILAC , INflA , IN60) 

IF  (MINJM  , QT i  1 NG A ) 60  TO  H33 
C  6€NCRAT£  CXTCRNALLT-LOADtO  NOOKS  ( N l 0 0 )  AND 
C  MAGNITUDES  OF  LOADS  ( CNLD ) 

IF  (NSS  ,EQ.  0 )NSS»NS/2 
N£«NZ* 1 

ND* ( NT* i  )  *nE*NR 
CLLD«-ftOO. 

IF  {NWOF  ,EQ«  1 )CLL0*-1000, 

CRES«0. 

NF*»NT 

if.  iHZ  ,EQ,  0)60  TO  IS  -  _ _ 

DO  1ft  m.NZ 
NP*NF*NT*l 
NLODm«NF*ND 

CDUM«.tt(ZC( m )-Ze(t3 ItCLLD 
CNLDI  I  )  ■CDUM*CRES 

CRCSaCDUM  _  _ 

18  CONTINUE 

CNLD(NC)«CRES 

14  IF  (NZ  ,EQ,  0 ) CNLD ( 1 ) *CLLD 
NL0D(NE)*ND4NF*NT*1 
IF  (NNOP  ,NE •  1)60  TO  22 

DO  20  IaI^JE  _ __ 

NL0D1 1*NC)«NL0D< 1 ) ♦NSS 
20  CNLDI  I«NE)rCNLO(I ) 

NE»NE«2 

22  NLOD ( NC*  1 ) «0 
NC1«NE-1 

....If  <NEJ_,Cfli  0)60  TO  2ft  _ 

00  24  Ul.NEl 
NLlPNLODlI) 

1 1 

00  2*  J«Ii |NC 
NLJ»NLOD< J) 

IF  I NOER  IN1.J)  ,6T.  NO£R  (Hi,  1)  |40  TO— tft- _ 

T 1 »NLOD I I ) 

T2SCNLDU) 

NLODd  )«NL0D(J) 

CNLD ( I ) vCNLO ( J) 

NL0D1 J  >  «T I 
CNLDI J)«T2 
NLlvNLODU) 

24  CONTINUE 

2ft  PRINT  200, (NLODI 1 ) i 1*1 ,NE) 


170 . 

171. 

172. 

173. 

1  7  «♦ . 
ITS. 

174. 

177. 

178. 
17*. 
1*0. 
1*1. 
182. 

183. 
188, 
188, 

184. 
1*7 1 
188, 
188, 
180. 
Ml. 
M2. 
M5. 
188, 
188. 
184. 

187, 

188. 
188j 
200. 
201, 
202, 
203, 
208, 
208, 

T04. 

207, 

208, 
208, 
210, 
211, 
212, 
213, 
218, 
218, 
*14, 

217. 

ir*.~ 

218, 
220. 
221, 
222, 

223, 
228, 
22ft, 

224, 


NAY»(^T*1  )  •  (  N Z ♦  l  1 

C  8ENERA7E  BOUNDARY  NODES  (NffSO)  AND  MAGNITUDES 

c  or  srecirico  displacements  «sroi 

If  (NNOr  ,eq.  1)80  TO  30 
MO* ( NT • 1  ) • ( NZ+ 1 ) 

00  24  1  u  1  (HO 
N8S0( 1 )•! 

24  srou,i)«o, 

MR»NR* ( NZ*  1  ) 

Mr«MO*MR 

NORbHO-NT 

00  27  !•!,**  _  _  _ 

NNSO(MO«! )*N0R*NT*1 
NOR*NNSO(MD*1 ) 

NNSO(Mr«l  )bNDR*NS 
SrO(MO*!,l 18200, 

2?  spdimf*; , n«20o, 

NfPO*M/*MR 
00  28  l«l,NSrO 

srou  ,2)*Q, 

28  SP0U,3)«200. 

ir  ( NOOT  ,EQ.  1  .08.  NO  1 M  ,E0«  2)80  TO  8ft 
00  81  j«i,Nsro 

I  *NRSp  ( J )  .  ..... 

Ul*(  1*1  1/NAY 
J2*N A Y • J | «  1 
U3*( 1 • J2 ) / ( NT« 1 ) 

IF  ( J j  , EQ i  0)5P0< U|3)*0. 

81  CONTINUE 

PO  83  U1  ,NTNN  _  _  _ 

Jl*( 1-1 )/NAY 
J2*N A  Y* j ) •  | 

J3«( I  - J2 ) / ( NT* 1 ) 

J8*1-J2-J3*<NT*1  ) 

IE  (J3  ,NE«  0)80  TO  83 

if  m  .t«u  o  toRj  jeq,  NHto_M_A3 _ _ 

ir  <JI  , EQ •  0)80  TO  *3 
NSrO*NSpO*l 
NNSO(NSrO)«l 
sroiNjrp,  n.aoo. 

SretNSPO. 2)*200. 
srotNSr^jji.o. 

83  CONTINUE 
80  TO  8ft 
30  NSPO*NE 

00  32  J.J  ,Nsro 
NN50( Jl-NLOOC J) 

sro (u. l ) *200. 
sro(j,2).o, 

32  SPO( J, 3)000. 

ir  ( N60T  .EQ*  I  .OR.  NO  1 M  (CQ,  2)60  TO  38 
00  82  U*1  iNSRO 
1 *NN SO ( J ) 

Jt*C 1-1  )/NAY 
J2*NAY*J| ♦) 

J3«  C 1 • J2  )  / ( NT* 1 ) 

tr  ( J3  , EQ •  0 1 SRO ( U . 3 ) *0 1 
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22  7. 
22*. 
229. 
230* 

231. 

232. 

233. 
23* . 
2  31 . 
23*. 
237. 
23*. 
239. 
2*0. 
2*1. 
2*2. 
2*3. 
2**, 
2*5, 
2*4. 
2*7, 
2*0. 
2*9. 
210. 
25|. 
252. 
213. 
25*. 
255. 
25*. 
267, 
250. 
259. 
2*0. 
2*1. 
2*2. 
2*3. 
2**, 
2*5. 
2**. 
2*7. 
2*0, 
2*9, 

270. 

271, 

272. 

273, 
27*, 
275, 
27*. 
277. 
270, 
279. 
200, 
201 . 
212, 
253, 


C 

c 


c 

c 


02  CONTINUE 

00  0*  !■( ,NTNN 
Ji«i  i-i  j/nay 
J2-NAY*JI*| 

J3»<  T»J2)/(NT*n 
J**I-J2«U3-<NT*1  ) 

If  <J3  ,NE.  0)60  TO  0* 

IT"  (JS  ,CQ.  NR  •  And.  J*  ,  to,  0)60  To  •* 

If  ( J  X  .50,  NR  *  AND ,  J*  .EQ,  NSS)60  TO  0* 

NSPD-NSROM 

NR9D ( NSRO ) *  1 

SRO (NSRD, I )  *200 , 

SPD(NSPD,2)a200. 

SPO  <  NSPO ,  3 ) *0, 

0*  CONTINUE 
05  NS!*NSPD"1 
00  29  I  a  I . NS  I 
NR  1 *NNSO  ( I  1 
l  1  •  f  ♦  I 


29 

3* 


DO  29  Jail, NSPO 
NR JaNRSD ( J ) 

If  (NOER(NRJ)  .67,  NOER(NR!))«0  TO  29 
Tl-NWSDU) 

T2.SP0U,n 
T3-SPD ( I ,2) 

T*aSPD  ( I , 3 ) 

NWSD( ! ) aNRSD ( J ) 

SPO (  !  ,  I )  *SPD  <  J » I  ) 

SPD ( 1 ,2)«SP0( J.2) 

SPO( I  ,3)-SP0( J.3) 

NRSD ( J ) »T 1 
5R0( J,J )«T2 
SPD< J,2)«T3 
SPD ( J  ,  3 ) aT* 

NW I aNRSD ( I ) 

CONTINUE 

RRINT  2  0  I  ,  ( NRSD ( l ) , l a | , NSpQ ) 

FORM  STIFFNESS  MATRICES  Or  TYPICAL  ELEMENTS 
STORE  IN  AUXILIARY  UNITS  I3,!*,t5 
IF  ( NCR  , EQ ,  1 ) 60  TO  *0 
MJNUMaNELEM/NER/3 

IF  (NELEM  ,6T»  MINUM-NER-3)MJNUM«M4NUM-1 
MJNAMaH|NUM 


AND 


JILL-0 

00  30  lal.NELEM.NER 
JILL«JILL*I 


IORUMa< JILL-1 )/MlNUM*i  3 

CALL  NQDESINOIM,) ,N00) 

NO  I aNOO ( | ) 

NQFrNOO ( NBa  I  ) 

CALL  INTE6(N0IM,l,NDl .NDF.SL) 

RR l TE  (l DRUM ) ( ( SL ( L 1 . L  2 ) • L I • 1 , NX ) »L2a I .NX) 

RRJNT  *00,I,(  (SLILI  ,L2)  ,Lla|  ,0)  ,1.2-1  ,5) 

30  CONTINUE 

*00  FORMAT  (25M0STIFF,  MATRIX  OF  ELEMENT , I */ ( 0E l I . 5 ) J 
BUILD  UP  6L0BAL  STIFFNESS  MATRIX  AND  SOLVE  FOR 

nodal  displacements 
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209, 

205 , 

206, 
207, 
200, 
209, 
290. 
29|. 
2  92, 

293. 
299, 

295. 

294, 
297. 

296. 

299. 

300. 

301. 

302. 

303. 
309, 

305. 

306. 

307. 
300. 

309. 

310. 

311. 

312. 

313. 
319. 

315, 

316, 

317, 
310, 
319, 
320* 

321, 

322, 

323, 
329, 

325, 

326, 

327, 
329. 

329. 

330. 

331. 
332* 
333. 
339, 
331, 

336, 

337, 
330, 
339, 
390. 


90  CALL  FACTOR 

00  9 |  l.l.NCLEM 

91  6RC  ( 1  I ■ 1 , 

00  80  M 1 ■ 1 ,NCTC 
REWIND  12 
M  X  •  0 

PRINT  202, Ml 

crlf«o. 

HINUM.NTNN 

00  93  l-I.NELEM 

IF  (BRC  m  «LE»  0, ISO  TO  93 

CALL  NODES  1  NO  I H , I ,N00 j 

N 0  I ■ NO 0 (  1  ) 

N0F«N00(NB>1 1 

call  coorojindi ,rqi ,tti .221 > 

CALI  CO0RO3(NOF,RDr  ,TTF  ,2zn 
RO* . 5 • ( RDF*RO 1 > 

R5".5»(RDF-R01 1 
2S*.5«(ZZr.2Zi ) 

IF  C TTF  ,LT.  TT1 1TTFP6. 283105307 
T0«,5«(TTF*TT1 1 
TS«.5*(TTr«TTl 1 
CALL  STRESS! 1 ,NOO,STRS) 

CALL  PSTRESCNOIM.LRM.STRS.R) 

IF  (LRM  , NE ,  |)S0  TO  31 

A  1 • ALF I »T0 

CRT  M , J ) «CS2*S 1 N ( A 1 ) 

CRT ( 2 , 1  J«.CS2«C0S(A1 ) 

CRT  C  1  ,  2  )  «C33*C0S ( A  1 >*SN3»Sn2«SIn(A1) 
CRT(2,2)«CS3*S)N(AI )»5N3*5n2*C0S(AI ) 

IF  CNDlM  .EQ*  2160  TO  1 1 

CRT  C  1  .3|pC33*SN2*S|NIAI )*SN3*C0S( A  I ) 

CRT(2,3)«-CS3»SN2*C0S( A| ).SN3«SIN(Al ) 

11  DO  23  I  A>]  ,  NO  1 H , N  J  D  A 
00  23  !  8* I , NO  1 M 
0UH»0 , 

DO  21  K«|,NQIH 

21  0UM«0UM*CR1F,IA)*CRT|K,18J 

C5N(  1  A,  1§)*DUH*0UH 
23  CONTINUE 

FOL  l»CSNf NO  IN, n*C2*C2*CSN(N01H,2 )«Cl*Cl 
F0L2«CSn(I ,1)»C2*C29C5NII ,21»CI»Cl 
IF  C  NO  1 N  ,CQ.  2)60  TO  12 
F0LI«F0l1»C3*C3*C3N(ND1M,J) *PAL*PAL 
F0L2«F0L2»C3*C3*C5Nf  j ,3>*PAL»PAL 

12  FOU  ,3)?Ral0/SQRT(R0L2) 

F0( I , I )«palo/sqrt<rol 1 ) 

CCCT»FO( 1  ,] )*TCRT 
CCCSbSCRT* ( FD ( l .31**0 f 1.1)) 

60  TO  33 

31  CCCT-FOI 1 .3J9TCRT 
CCCS«2.»F0( I ,3)»SCRT 
33  RT|«*P(  i  )/FOU  ,3) 

RT2«F(N01H)/CCCT 

RT  3« ( F ( NO t  M ) «P ( 1 ) l/CCCS 

BRCC I IBAHAXI IRTJ .RT2.RT3) 

IF  IBRC(I)  ,6T,  CRLF ) CRLFbBRC ( I  ) 
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3*1, 

342, 

343. 

344, 

345. 
344. 

347, 

348. 

387, 

380. 

381. 

352 , 

353. 
364, 

355. 

356. 
317. 

388, 

357. 
340. 
34|. 

342. 

343. 

344. 

345. 
344  • 

347. 

348, 
344, 

370, 

371, 

372, 

373, 

374, 

375, 
374, 

377. 

378. 
374, 
380, 
38], 

382, 

383, 

384, 

385, 
384, 

387. 

388, 
384, 
340 , 
341  , 

342, 

343, 
34*. 
345. 

344, 
347, 


43  CONTINUE 

C  DETC8HINE  ELEMENTS  THAV  ,'aILED,  LA  (  I  } 

DO  44  l-l.NELEM 
I F  < QRC (!)  ,LC.  0, ) GO  TO  44 

IF  (BRC(I)/CRLF  *LT»  (1«*TRN))G0  TO  44 
I  1  •  I  1  ♦  1 

L  A ( 11)4! 

BRC  ( !)*•]. 

44  CONTINUE 

C  CRITICAL  LO AO 
CRL» I  000 ,/CRLF 
PRINT  144.CRL 
N 1*NTNN*N0ER(NS'  1  )M 
N2«NTNN.N0ER(NST2)*1 
M1«NTNN*N0ER(MST1 )♦{ 

M2»NTNN«N0ER(M5T2)*1 
TEMP*(RC(2)-RC( 1 ) )*CRLF 
STRN* ( 0 (  I  ,N2 )-0( I ,N | ) ) /TEMF 
STRm*<0<1 ,m2)«0( 1 , M 1 ) ) / TEMP 
AVE*,5«(STRN+STRM) 

MR  I  NT  203 
PRINT  204.STRN 
PRINT  205.STRM 
PRINT  204, AVC 
RERIND  13 
REWINO  14 
RERINO  15 
NH( 1 )«C 
NM(2)*0 
NM( 3)*0 
PRINT  207 
NLF*(NE*1 )/2 
DO  42  1*1 ,NLF 
NLO*NLOD ( 1 ) 

NUN*NTNN»NOER(NL0)7l 
CNL1»CNL0U  MCRL/CLLO 
021*0(1 , NUN ) *CRL/ 1 000 , 

IP  (NE  ,EQ,  l ) GO  TO  45 
J*I*NLP 

IP  (J  • GT ,  NE ) GO  TO  46 
NLA*NL0D( J) 

NAN«NTNN*N0ER(NLA)«1 
CNL2*CNL0(J)«CRL/CLLD 
DZ2*D (  1  ,  NAN )»CRL/ 1000, 

PRINT  204,NL0,CNL1  ,DZ1  ,NL A. CNL 2,1)22 
GO  TO  42 

45  PRINT  204.NL0.CNL1 ,0Z1 
42  CONTINUE 
PRINT  211 

IP  (NDlM  • EQ •  2  )  PR  I  NT  214 
IP  (NDIM  , Efi .  3  )  PR l NT  212 

C  AN°  ™CN  DC0UCT  5tIPPNE5S  OP  FA1LE0  ELEMENTS 

DO  74  1*1  ,  i  1 

c  shear  failure 

I C0DE*3 
LDA*L A  (  i  ) 
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398. 
39*  , 
•♦0(5, 
“01, 

902, 

903, 
909, 

905, 

906, 

907, 

908, 

909, 

910, 

911, 

912, 

913, 
919. 

915. 

916. 

917. 

918. 

919. 

920. 

921. 

922. 

923. 
929. 

925, 

926. 

927. 

928, 

929, 

930. 

931, 

932. 

933. 
939, 
935. 

934, 

937, 

938, 

939, 
99n. 
991  , 

992, 

993, 
999. 
995. 

994, 

997, 

998, 

999, 
950. 

981. 

982. 

983. 
989. 


C 

c 


CALL  N00ES  (ND1M.L0A ,NOD  ) 

NQ I «NOD ( 1 ) 

NOF«NOD(NB»| ) 

CALL  C00R031N01 ,ROI ,TT!  ,Zzi  ) 
CALL  C00RD3(NDF,RDF,TTF,ZZF) 
R0«.5*(R0r*R0I I 

R5*.5* (ROr.ROl ) 
ZS*.S«(ZZF»ZZ1  ) 


IF  (TTf  ,1.7,  TT  |  )  TTF«6  •  283  1  85307 
T0«.5«(TTr*TT! I 
T5«,8*(TTr.TTl ) 

CALL  ST8ESS<LDA,N0D,STR5 > 

CALL  FSTRE5(N01M,I ,STRS,P) 

IF  (LRM  ,NE.  1 (GO  TO  35 
FT  <2)«Fd(L0A, 1 )*TCRT 


2i"l**<PC><LOA#3,,fcPD<LOA.l  )  )»SCRT 
GO  TO  37 


35  PT(2)«RD(LDAi3)«TCRT 
FT (3I«Ro(LOA,3)«5cRT 
37  FT ( 1 ) «Fo ( LD A , 3  I 

FT1«-R(  |  )  /PD ( L  0  A ,3 ) 

RT2«P(N01M)/RT(2) 

rT3«,5*(P(nD1M)»F(1))/ptI3) 

IF  (RTl  ,GT »  RT2)G0  TO  96 
TENSILE  FAILURE 

IF  ( RT2  , GT ,  RT3 ) I  C0DE*2 
GO  TO  97 


COMPRESSIVE  FAILURE 
96  IF  (RTl  , ST ,  RTS ) t CODE* 1 

.  <mu  ,  .'MniiDe*; Um'  '  >L,,‘ >N' 

IF  (MI  , EQ ,  NCYOSO  TO  80 
IF  (HER  , EQ,  1 (GO  TO  80 
NN*(L0A«1 ) /NER* 1 
nmn»(nn.i >  /m  i  nam* i 
IORum»nmn*i 2 

NAM.NN-(NMN-1 (•MINAM-NM(NMN) 

IF  (NAM  , EQ ,  0 ) GO  TO  52 
00  98  IR«1,NAM 


NM(NMN(.NN.(NMN-1  )  *'  M  I  n  a  m 
GO  TO  62 

50  CALL  INTEG(N01M,L0A,N0I ,NfiF 
52  IF  ( 1C00E  , NE ,  2JG0  TO  60 

IF  (NTEP  ,EQ.  1 ) Go  TO  40 
NW«1 


,SL) 


IF  (ISO  ,EQ,  1 (GO  TO  54 
1  SO*  1 


e l *FD ( LO A , 1 ) 

C2«E  1 

IF  ( NO  I M  ,EQ«  2 ) E2*0« 
C3«0« 

F1*P0(L0A,2> 

F2«P1 


i 


F3*P1 

G 1 «E 1 /2 , / ( l . «p | ) 


64 


*35. 

*<34, 

437. 

438, 
454, 
440. 
44|  , 

442, 

443, 
44*4  , 
448 . 

444, 

447, 

448, 
444, 
470. 
47  I . 
4*2. 
4  f  3 . 
474, 
'+75. 
474. 

477. 

478, 
474. 

480. 

481. 

482. 

483. 

484. 

485. 

486. 

487. 

488. 
484, 
4»0. 
4»l, 

442. 

443. 

444. 
448. 
444. 

447, 

448, 
444, 
S00. 
tOl, 
902. 

803, 

804. 
808, 
804. 

307, 

308. 
804, 
810* 
«U. 


C2-G  » 

G3«G  1 
GO  TO  38 

c  for  originally  anisotropic  discs 

54  CONTINUE 

c 

58  CALL  1NTEG(ND!M,LDA,NDI .NDF.SLD) 

DO  34  L  1  • 1 , N X 
DO  54  L2a 1 f  NX 

54  SLDfLl *  L  2  )  >5LD ( L l ,L2)-SL(L1.L2> 

GO  TO  42 

40  DO  4!  L  1  •  1  ,NX 
DO  4  t  L2*  I .NX 

*1  SLDIL1 ,L2)*-40SM»SL(L1 ,l2) 

42  DO  77  L |*|  ,NB 
L I L*NOD ( L  1  1 
LX*N0CR(L1L>«MX 
MX*N0ER(L1L) 

IF  (NX  ,LT.  MlNUMjMINUMaHX 
IF  (LX  ,LC.  0)G0  TO  45 
IF  (LX  ,EQ.  I  ) GO  TO  48 

lx-lx-i 

DO  44  I  A* I  .LX 

Read  <  i?. inn3,((S(k,jjik»i.ndim},j*i,nn3) 

44  CONTINUE 
GO  TO  48 

45  NXXal-LX 

DO  44  I  A* 1 .NXX 

backspace  ji 

46  continue 

48  READ  < 12)NN3,< (5(K,U) ,K*J  ,ND1M)  ,j«i  ,NN3) 
BACKSPACE  12 
I DUM* ( L I •  1  )  *ND  I  M 
DO  44  I  A* 1 .NDIM 
I  I  An  I  A* I DUM 
DO  44  IB»WND|H 
I  I B • 1 B ♦ | DUM 

44  S(IA,IBIpS(IA|1B)«SLD(I1A,|IB) 

IF  (MX  . EQ«  NTNN ) GO  TO  74 
NN*NN3/N0IH«1 
00  75  L2* 1 . NB 
LAL4N0D(L1) 

JXaNOCR ( L AL ) 

IF  (JX  ,LE.  MXIGO  TO  71 
NiCaNBRC (MX )•! 

DO  71  L3«l ,NN 
LL3«L3 

IF  ( LOC ( NBC+L3 1  .CO,  JXJGO  TO  72 

7 1  CONTINUE 

72  I  DEM«LL3»ND 1 M 

I D AMa ( L2» I ) • N  D I M 
DO  79  1  A  ■  1  , ND I M 
I  I  A ■ I  A  *  I  DUM 
DO  73  !Bal,NDlM 
1 JAa  t  B* 1  DEM 
! !■•! Ba I  DAM 

7  3  S< 1  A  a lJA>aS( IA|I JA)«9LD( I  I  A, I  IB) 
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512. 

*!3, 

SIM. 

SIS. 

su. 

S  1  7 . 
5l«, 

sit. 

S20 . 

52  1  , 
522. 
52). 
52m  , 
525, 
526  • 
»27. 
5  25. 
S2t , 
530. 

53  i  , 

532. 

533. 
5  3m. 

535. 

536. 

537. 
S3  A , 
53t, 
SMO. 
SMJ  . 
SM  2, 
SM3. 
5MM, 
5MS. 
5M6, 
5M7, 

5  M  A , 
SMt. 

550. 

551, 
552  , 
553, 
55M, 
SBS . 
566. 

557, 

558. 
55f. 
560, 
56|. 
562. 
56  3. 
56M  , 
665. 

566. 

567. 

568. 


75 

76 


77 

7 1 

80 

100 

102 

10M 

106 

ltt 


» J"l , NN3) 


.NN3> 


CONTINUE 

. . 

(  12)NN3,( (S(K,J| ,K*| ,N01H) ,j«t , 

CONTINUE  *' 

l  SO-LRM 
conti nuc 
call  fc  tor 
CONTINUE 
FORMAT  (9)5) 

format  drjo.oi 
FORMAT  (MF10. 0,215) 

. w.  format  inoi 
ltt  FORMAT  (17H0CR1TICAL  LOAO  .  .E12.A1 

200  FORMAT  ( 27H0EXTERNALLY-L0ADED  NODES 

201  FORMAT  (J8H050UNDARY  NODES  .. , ) 0( 1 M  iH.i 

2/(18X,IO(IM,lH,m  •10(1“, 1H») 

202  FORMAT  (IOhOCYCLE  N0.,13»5h  •  •••) 

»i . 

“  « !swr  :  i 

.  NODE', 5*. ♦  LOAD  •,5x,»01SRLACCHCNT',I0*)) 

r2SM*I  (2{  ,‘,»*,*'E*0'5*6X,Cl0.5,10in 

211  FORMAT  UHO.tX.ltHSTRESSES  IN  FAILED  , 

I  3 ] HELEMENTS  (DUE  TO  1000-LB.  LOAD)) 

213  FORMAT  ( I  6 , MX , A6 , A5 , 1 0 ( | X , e I  0 • 5  )  ) 

212  format  element  fail.  mode  ra 

i,  fXlAL  SHRRC  SHRRZ 

2'  F  R  N  c  R  PRNCP2  PRNCP3 

2|M  format  (.  element  Fail,  mode 
I  '  SHRRC  PRNCP 1  PRNCP2 

FORMAT  (20HIMATERIAL  PROPERTIES) 

FORMAT  ( 23m  ELASTIC  MODULUS  .  ,|PE|2,6) 

foissons  ratio  ■  ,fm,2) 

ALLOW,  COMPR,  STRESS  ■  .2PE10.3) 

*  I  PE  1 2*6  , 


6 T  LOADED  NOOES*/ 


300 
30  1 


RADIAL  CIRCUM, 

SHRCZ 
ALLOW. • ) 

RADIAL  CIRCUM. 

ALLOW, * ) 


302 

303 


306 


(22H 
(  28 H 

(  23H 


FORMAT 

format 

FORMAT 

I^H  TO  ,  I  PE  1 2 , 6 ) 

307  FORMAT  (22 H 

|NH  TO  ,  FM , 2 ) 

FORMAT  (28H 
|RH  TO  .2PEI0.3) 

format  (i»H  elastic  modul i /ex ,shei 

18X.5HE2  -  ,IPE|2.6/8X,5HE3  -  . 1  PE  1 2. 6 ) 
FORMAT  (20H  POISSONS  R A T I  OS /8 X , 5hP I 
1  X.5HP2  ■  ,  FM , 2/8x , 5hP)  ■  ,FH«2) 

shear  modul i /ex ,shgi  ■ 

uj,;1 


308 

309 
3|0 
31  I 


elastic  modulus 
poissons  ratio 
allow,  compr . 


*  PM , 2 , 
STRESS  ■ 


, 2PE I  0 , 3  , 

■  , I  PE  I  2 , 6  / 
■  , FM  »  2/ 

, IPE12.6/ 


I17H  of  ALLOW.  roMMB.i  '  n»«» 


305 

313 


I  I  7H  of  allow,  compr, ) 
FORMAT  ( 28H  ALLOW, 

1  I  7H  of  allow,  COMPR.  , 
FORMAT  ( 28hOD I  MENS  I ONS 


TENSL,  STRESS 
SHEAR,  STRESS  ■  , 
OF  DISC  SPECIMEN) 


pm  ,  2 , 

FM.2, 


3|*  FORMAT  (UH  DIAMETER  ■  ,F5,2> 

315  format  < i 7 h  thickness  .  ,fs.2> 

316  format  UtMORADlAt  COORO  l  N  ATES  /  (  *♦  X  » |  0  ( F6 . 2  ,  I  H  ,  )  )  ) 

3|?  FORMAT  <38M0CIRCUMFERENT J al  COORDINATES  1  DEGREES ) 

l/(RXilO(F6t2ilH,))) 

318  FORMAT  (  1  RHOZ'COORO 1  NATES/ { NX , I  0 1 F6 «2 • 1 H t )) ) 

31*  FORMAT  ( 2RM0T0T Al  NO.  OF  ELEMENTS  •,!•*> 

320  FORMAT  ( 2 1 HOTOT  AL  NO,  OF  NODES 

326  FORMAT  <♦  AN J SOTROF  l  C ♦ /8 X , *  ALPHA  1  ■  ',2PEI0,3, 

1*  DEQREES*/8Xi,ALFHA2  •  •  »2PE10.3/8X,*AlFHA3  •  *»2pE10»3> 
R33  STOP 
END 

END  OF  COMPILATION!  NO  DIAGNOSTICS, 
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56  9 « 
570, 
57|. 

572, 

573, 
S7<4. 

575, 

576, 

577, 

578, 

579, 
510, 
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7  • 
3  • 

f  . 
S  . 
A  • 
7. 

n  . 
o . 
I  n . 
I  1  . 
I  7 . 
I  t . 
l  r  . 
I  s . 

1  A  . 

I  ■' , 
1  -I . 

1  9. 
7 (I . 

2  !  • 
77. 
?  ) - 
2'i . 
?G  . 
2A. 
?7. 
7.  n . 
29. 
3n, 
3| . 
37  . 
3  3  . 
3R  . 

is* 
3  A  • 
37. 
3fl. 

3  V  . 
so. 

*♦  l  • 
*(;•. 
*n. 

M'l . 

hs. 

*<a. 

H  7  . 

hr  . 

H  9  . 
Sn . 
Si  . 
S7. 
S3. 
SH. 

ss . 
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GENERA TrS  NOOKS  C  n  N  NF  C  3  F  (>  TO  NOOK  I  AND  ARITES  THEM 
IN  AUXILIARY  STORAGE 
KOR'IS  KMfirH  COUNfFR  ARRAY  ,  NRRdll 


rOMMOrf/R.i  r.l/LOr  (  2snoo  )  .NTflN  ,N0  ,NR  ,m  STY 
rn.lMON/rNNKV/  IRRC  (  SOP  |  .nHFO(SOO)  .NOFRISOO) 

COMMON/  iflFrlf/NR.NS.N/.NWOp,  ,r>|M 
0  1  mens  MM  MCON(.?/,  ,|.s.M  TNr,A  ,  INGO) 

N  T  a  NS  -||  ,f>3P 

fIO«  <  NT*  )  )  •  (  HZ  ♦  1  1 

M  =  -  l 

00  1A  t  a  1  ,  »j  T  N  N 

M  =  M  ♦  | 

1  1  =  <  I  “  1  WHO 
I  ? 3  MO • J  I  *  I 
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'ICON  <  3  MNCON  (  ?  I  ♦  1 
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N  c  0  M  (  6  )  a  J  ♦  M  o 
I J  C  0  N  (  0 1 *NCON  (  1  ) ♦no 
N  C  0 N  ( 7  |  aNCON ( H ) ♦  NO 

IE  I  1 1)  1 M  «F'i*  7)5o  TO  3 
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'ICON  (ll]|  =NCOf|  (  7  )  -NT-  1 
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0  C  0  M  (  1  2  J  »  N  C  0  N  (  1  I  )  ♦  N I) 
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ICON  (  I  0  )  a  NCON  (  1  I  ♦NT  ♦  | 
ocon ( i9i»ncon(7)*nt+i 
'J  C  0  N  (  2  0)  a  N  C  0  N  (  3  )  +  N  T  ♦  I 
'ICON  (211  aNCON  (  ? Q  )  *140 
'ICON  (  22)»NC0M(21  I  ♦NO 
■  J C  0 N  (  2  3 ) a N C ON ( A ) ♦ N  T ♦  1 
ICON  <  25 ) *NC ON (  I R I ♦NO 
He  ON  (  2S  )  =NCON  (  25  1  ♦Nil 

■‘I  C  0  N  (  2  6  )  a  J  ♦  N  T  ♦  1 

II  M3  •  £Q .  0  .OR.  J  3  .  F  Q  .  <g  Z  )  C.O  To  1 
GO  TO  3 

I  UI*<J3/nZM)*? 

N  F  a  N  I  ♦  p 
00  2  J  a  N  I  ,  N  F 
?  'ICON  (  J  I  a  l) 

1  11  ( J  1  •£-').  o  .  0»,  j|  ,F'i.  nr  )  f,0  Tp  H 


S  A  • 

(in  Tn  a 

S7  . 

9 

N  l  =  1  *J l /  N  rt  «  9 

S  P  . 

N  r  *  N  l  «■  2 

5  9. 

00  r>  jsfj  |  ,  N F 

ao . 

N  C  0  N  (  J  )  a  0 

a  i  . 

N  C  o  N  <  J  +  a  >  *  0 

a? . 

5 

NC0N(JM7)*n 

M. 

A 

tr  (N;jof>  .pr).  i  i g o  to  ii 

A  9  . 

IF"  (  M  .  r  rj  .  n  )  (,  J  To  7 

• 

I  r  t  M  .  r  :J ,  NT)  & 0  TO  g 

A  A  . 

GO  TO  11 

A  7  . 

7 

NC  ON  (  1  i  30 

ah. 

NCOM  (11)  *0 

A  9  . 

N  C  0  N  (7)'!) 

70. 

N  C  0  N  (9)s0 

7  I  . 

N  C  0  N  (  t  A  )  a  0 

77. 

NCON  (  1  !i  )  »IJ 

73. 

J  C  0  N  (  I  A  )  «  n 

79  . 

NCO M ( 2 ft  >  * n 

7c,  . 

•JC0N(79)3U 

7  A  . 

so  TO  11 

77. 

n 

Of)  9  J»i,g 

79  . 

N  C  0  N  (  J  )  a  0 

7  9  . 

N  C  0  N  (  J  ♦  9  )  a  n 

no . 

9 

NCON  ( J ♦ |  7 ) a  ft 

»  1  . 

M  a  -  1 

7  . 

)  t 

NI39C  (  I  )  »?A 

«  3  . 

IK  (  1  0  I  M  .  EG)  .  2  )  Ng9c  (  1  )  a8 

AS  . 

0  M  a  N  0  K  C  (  1  ) 

«S  . 

NCON <  MM+  l ) a0 

A  A  . 

M  M  1  a  *1 M  -  1 

M  7  • 

0  0  |9  J a  1  , N  M  t 

nn . 

IK  (NCON(J)  •  NFC  •  0)  C,0  TO 

A9. 

>)  l  *  J  ♦  I 

9n, 

no  17  K  a J 1  ,  MM 

9  I  . 

K  Ja< 

97  . 

ir  I  N  C  0  N  (  K  )  «  NF  «  0)  GO  TO 

93. 

l  7 

C  ON  T  1  Ntlf 

9  9  . 

N09C ( I ) a J- J 

9C,  . 

G  0  TO  is 

9a. 

t  3 

NCON (  J ) a NCON ( K J ) 

97  . 

N  C.  0  N  (  K  J  »  a  0 

9  n. 

l  9 

CONT 1 N jr 

99. 

1  s 

N  •)  sMKHC  (  I  ) 

)  fio . 

0  0  in  J  a  i  ,  'j  n 

1  0  I  . 

I  0 

1.9  0  (  J  ,  1  )  a  »CON  I  J  ) 

107. 

l  A 

CONT I  NOT 

103. 

H  i;  r  1 1  s-  tg 

ton. 

F.  N  0 

rNn  oi 
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i  . 
?  • 
3 . 

‘4  . 
5  . 
A  . 
7  . 
n . 

9 . 

I  n . 
1  I  . 
I  ?. 
l  o . 
1  ** . 
I  s . 

1  A. 
1  7  . 
1  1 . 
I  9  . 

2n  . 
?  t  . 
??. 

?3. 
2ft. 
2ft . 


ftP  . 
ft9. 
5r . 
5-1  , 
52. 
50. 
5*4  . 
55. 


C 

c 

r 

c 


COMMON /rcurn/LOe  (  25000)  .nTnn.NIJ.NP.MtST  V 

C0M,'4<vi/rUNKY/NnRC(Son),f1tRE-n(sn0),NOF,<<5nn) 

COMMON/  )NF  F  R  /  Nl<  .  N  5  .  N  Z  iNWOF*,  J  ft  j  m 
!)  |  Mr  MS  1  ON  L  *  1  1  0(1  )  ,  LSO  (  I  NGA  ,  I  NGO  ) 

KE'VIMO  |t 
print  ion 
j  =  n 

0  Z  E  T  *0 
no'iMn  i  non 
00  7  ) = 1 , NTNN 

I  F  (  MORf  (  I  )  .  |.  T  .  n  )  00  TO  7 
IF  (NnNCUl  •  GE  •  NDUM)  GO  TO  7 
NOUMbNHRC (  J ) 

•J  L  0  C  T  =  l 
COMT 1 NOF 
•JLOC*  Jl.OCT 
J  L  0  C  T  =  1 1 
J  a  J  •»  1 

NR  E  0  ( .1  )  a.ILOC 
OOFR  (  Jl.nC  1  n  J 


2  A  . 

M  N  ■  N  ft  N  C  (  ILOC  I 

27  . 

'VR  l  T  ft  (  I  1  ,  /no  I  NN,  (LSO(L  ,  JLOC  >  .1  ■  1 

,NM  ) 

2ft. 

N  7  F  T  a  N  Z  F!  T  *  N  N 

29  . 

NMsH  N 

30. 

II  (J  .EO.  NT NO- 1)60  TO  2ft 

3  1  . 

Of-  ?  3  1  a  1  ,  NM 

3?. 

ELOOLSO  (  1  ,  JL.OC  I 

3  3  . 

4  MaMHfsC  ;  K1.0C  ) 

3*4  . 

1  r  (MM  . G  T ,  INGA)  GO  TO  3U 

3  5 . 

0  0  !  *4  J  J  c  1  ,  M  M 

3  A  • 

J  M  r  J  J 

37  . 

ii  nsnuj.KLOci  .to.  jLOcir.o  Tn 

15 

3ft. 

1  ft 

CONTI nur 

39  . 

1  5 

M 1  a  MM  -  1 

on . 

(  0  l  A  l.l.  s  Jfl  ,  M  1  ,  1 

ft  l . 

1  A 

i  son  i.  >  k  t.  o  r  j  *  i.  s  o  ( u  i  ♦  j  ,  k  l  o  c  i 

ft/ . 

MftWf  (  KL  OC  )  «  N  ft  R  C  (  K  L  0  C  )  -  i 

‘13. 

M  M  a  N  ft  P  C  (  K  L  n  r  I 

ftft  . 

0  0  19  Jv|a  I  ,  MM 

'<5. 

IF  (  LSD  (  J  J  ,  jL  f*c  1  .  EQ.  KL0C)6O  tq 

1  9 

ft  A  . 

On  1  H  L  L.  a  1  .MM 

ft  7  . 

U  1L5|)(.IJ,  JLOC)  •  FO ,  LSD(U  ,KL0C)  IGO  TO 
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NftRC  (K.Loe  )  tMHRC  <  iq  OC  )  ♦  I 
M  =  NORC  <  K  L  0 C  I 

LSO(M,KL  OC  )  aLOOl  JJ.  JI.OC  I 
19  FONT!  MIC 

ir  (NHRC(KLOC)  »  G  T  .  NM)  GO 
NM  =  MlRC  I  K'LOC  ) 

JLOrTaKl.  OC 


I  9 
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(  s . 
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7n. 
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7  A 

conti nor 

S‘i  . 
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l.  AKsl  A  (  r  1 

P7  . 

NS  J  a  US  J  *•  1 

fin . 

?  7 

i or i  vi  j ) » nor R <  l \k ) 

A*?  . 

NS./aMS  J*  I 

9n . 

L  0  C  (  N  S  J  »  a  f 

•>  1  • 

7  9 

r  (i nT  J  N 0 F 

9?  . 

PRINT  1  0.?  ,  MTNH  ,  NRFO  (IJTNN  » 

9  3  . 

PRINT  | 09 , PRC  r 

9 '4  . 

ion 

‘ORMAT  (*1  RESULTS  OF  nOOAL  REORDERING 

9S  . 

! 

•  o  n  f  w  n  o  •  n  R  i  r,  •  no.  nodal  connections 

°  A  • 

ini 

format  (•o**«»*  array  length  is  excffofd  * 

9  7  . 

20'1 

FORMAT  {  .3  2  I  ‘4  ) 

9  n . 

in? 

FORMAT  (T6,PX,ta,<4X,221S/(2|X,22IS>) 

99  , 

1  0  '4 

1  H'iiaT  (  3  ?  HOP  E  RC  Em  T  A  GE  OF  nON-ZFRO  TERMS  a 

inn. 

GO  TO  3  | 

I'M. 

30 

i  I  ST YsMrt 

1  nr . 
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RF  TURN 
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End 
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1  . 
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3. 

H  • 
5. 
6  • 
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13. 
19. 

15. 

16. 
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SUBROUTINE  ELTOnO  (  10  I M  ,  I ,LCOn ) 


C  GtNEHATES  ELEMENTS  CONNECTED  TO  NODE  I 

C 

c 

COMMON/ JNFER /NR ,NS,NZ ,N*OP,NDlM 
DIMENSION  LCON(B) 

NT  *NS-NrtOp 

NO* <  NT*  1  1 • <  NZ*  1  ) 

NG*NS*NZ 

IF  (1DIM  . Eg .  2)  N  G*  NS 

J  1  ■ ( 1  -  1  ) /ND 

J2«N0* J 1 ♦ j 

J3* ( l-J2»/<NT*l ) 

K  J*  J3 

1 F  1 NWOP  .Eg .  1)  K3*0 
LCON  <  l  )*1-(nD-N(,)*j1-K3 
LCON ( 2 1 «LCOn  <  1 ) -  1 

IF  1NWUP  .Eg.  I  .AND.  1  •  EQ  •  ( J2 ♦ J3 *NS »  » L C ON ( 2 ) «LC 0N 1 1 » ♦ N 1 
LC0N(3)*LC0N(2)"N(, 

LCON  <  H ) *LC  ON (  J  )  -Ng 
IF  (  IOIM  .Eg.  2 ) Go  TO  3 
LCON (5J-LC0N ( 1 )-NS 
LCON ( 6  5  *LCON ( 2 1 -NS 
LCON ( 7 )*LCON ( 6) -Ng 
LC0N(8)*LC0N<5)-Ng 

IF  <J3  .Eg.  0  ,0H,  J3  ,EW.  NZ)  GO  TO  1 
GO  TO  3 

1  Nl*5-J3/N2*s 
Np *N  I  ♦ 3 

DO  2  J*N I  ,Nf 

2  LCON ( J 1 *0 

3  IF  <01  .Eg.  0  .ON.  Jl  *E0.  NK>  GO  TO  9 
GO  TO  6 

9  NI*3-JI/NR*2 

N  F  c  N  I  *  I 

00  5  0 ■ N  J  ,NF 

L  C  0  N  (  J  )  ■  0 

5  LCON(J*H)»0 

6  I F  ( NWOP  .Eg •  I  ) GO  TO  9 
J  9 ■ 1  -  J2- J3* ( NT* 1 ) 

IF  <09  .Eg.  0)  GO  TU  7 
IF  ( Jn  • E  Q .  N  T 1  Go  TO  8 
GO  TO  9 

7  LCON ( 2 ) *0 
LCON ( 3 1 *0 
LCtN ( 6 ) *0 
LCON ( 7 ) *0 
GO  TO  9 

8  LCON ( I )*0 
LCON  (  <4 )  »U 
LCON ( 5 ) *0 
LCON ( 8 1 *0 

9  IF  <  I  D 1 M  .Eg. 


2  I Gg  TO  IQ 


I 

I 
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S6. 

J  1  ■  L  C  0  N  (  M  ) 

S7  . 

LCONtH»-LCON<S) 

SB. 

LCON ( S ) ■ J l 

sv. 

J i «L  CON ( 0 ) 

60. 

LCON  <  3 ) -LCON ( 6 ) 

6  |  . 

LCON ( 6 ) ■ J 1 

62. 

10  REIOKN 

6  3. 

t  NO 

ln0  of  compilation:  no  diagnostics. 
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I 


l  . 
2 . 
J . 
'( • 
b  • 
6  . 
7  . 
b  • 
V  . 
1  L  • 
1  1  . 

1  2  . 
1  3  . 

1  *i  . 

1  b  . 
1  6  • 
1  7  . 
I  b. 

1  v. 
iU. 

2  1  . 
22  . 
?  J  . 
2M  . 
2b. 
2  6. 
27  , 

?e. 


c 

r 

C. 


SUfcKOiniUL  I'iUotb  (lulH.l.NOo) 


OLl.LKATtb  NUDaE  NuMbERS  OF  ELEKEnI  1 


COMMON/ J.IF  th/NK  ,  NS  ,MZ  *N*Op  ,  MU  1  h 
DIMENSION  NOIM8  1 
N  I  *  N  S 

If  lN»Of  ,t(j.  l  )  M*hS-| 

Mb*-  l  NT*  1  I  •  (  N  2  ♦  1  ) 

N  (,  t  N  S  *  N  / 

If  i  i  c  l  m  .  t  o .  2  i  f«  c.  - 1.  s 

i ;  - 1 » - 1 > /mg 

I  2  »  N  U  •  1  J  ♦  | 
l-JM  l  “ I  2  > /Nb 

K  J*  1  J 

IF  t  MO  0  F’  •  E  W  •  IF  K  J  “  U 
M  0 1  I  l  )  «  |  *  (  M»  -  M  o  )  •  1  l  ♦  K  3 
NOD ( 2 1 *NUO ( 1 1 « i 

IF  (NIF'OP  .Eo.  1  *  A  NO  *  I  .tO.  |  /  N  S  *  N  S  )  N  0  U  (  2  l  *  N  U  0  (  1  )  -  N  I 

NoOIJ»»NOOl2F«-Nu 

NO  l  (  b  1  «  m  0  D  (  j  )  ♦  N  U 

If  (  I  DIM  • E  U  •  2)00  TU  lo 

N  0 1  I  b  )  ■  1. 0  C  |  1  )  ♦  N  1  ♦  i 

N 00(61* NuO(2)4Nf4l 

N 0  1  (  7  )  »M00  <  t,  1  ♦MO 

MOt (6)*m0D(s1>N0 

ID  N  E  i  u  K  N 
EkO 


{.mo  of  comhilatjon: 


no  diagnostics* 
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1 , 

SUBROUTINE  COORD3I I .ROC.TTC.ZZC) 

2. 

C 

3, 

c 

assigns  radial, circumferential,  and  z  ordinates  to  noof 

• 

c 

S, 

COMMON /JNFER /NR, NS, NZ, NWOP.NO I M 

^  • 

COMMON/lEVTF/RC(3O},TC(S0>,ZCC3O) 

7  , 

zzc-o. 

fl. 

nt»ns-nwop 

9, 

NO* ( N  T  *  l  ) • ( NZ  ♦  1  J 

in. 

JI" ( I  -  1 ) /NO 

I  i  , 

ROC *RC ( J 1 ♦  I  ) 

12, 

J2*N0* J l ♦ l 

1  3. 

J3*  <  1  -J2 ) / ( N T ♦  1  ) 

1  9  . 

IF  ( NO  I M  , F 0 •  2 ) GO  TO  I 

IS. 

ZZC*ZC( J3*1  I 

1  • 

1  J9* I - J2- J3* ( NT*  I  ) 

1  7  . 

TTC-TC  <  J M ♦  I  ) 

1  A. 

RETURN 

19. 

END 

fnd  nr  compilation:  no  diagnostics. 


1.  SUBROUTINE  UN  I FRM | A , B , C > 

2.  O-RANUN(R) 

3  •  C«A*(B»AI*0 

m.  return 

S.  END 


fnc  of  compilation: 


NO  diagnostics 


bUttKUUl  iNt  hlKXbl  I(<lM.I>K»K,TZ,Z,t>b) 


^.OKhb  THE  MHKU  b  AND  rHE  MAfKlX  PRODUCT  0*b 


C jmhuN/sANNA/PD < HpO ,3)  ,FNCTR(b)  ,FNCTT  (fi)  ,FNCT  Z(a> 
C0MN0U/B(,ANL(/C(B»8),b(6,2m,H 
CohflON/7  ANA  1  /  hAD  »Rb  ,T0,  TS  ,ZS*NA  ,HNB 
CuMMUN/TATAY/ 1SO.A1  ,A2,AJ,E1 ,E2,t3,PI  ,P2*PJ, 

1^1  i 02  »  (i  3 

COMMON  / NaNaY  /  T  T  (  3,3)  •  N  Hi 

DIMENSION  »U<6.2‘<)  (FNCT(tt)  ,otM6.0)  l0(6,6>  ,T(6,6) 
Na  I  SNA- 1 

NY=4*(1DIM-1»«Iu1m 

It-  (1U1M  .  E  U  •  2)00  f  0  3 

»•  fJC  I  (  1  )  *  I  .  /K  Ao 

)NU(2)«K/kaO 

F  NC  T ( 3 )*T  Z/R AD 

F  r,C  1  (  <UsF  NC  I  (  3  )  «K 

F  NC  T  t  5  I  at  Z  /R/iU 

FnCT<6)*FncT(5)«N 

F  N  C  1  (  7  )  =  V  N  C  T  I  b  I  •  T  i 
F  NC 1  ( B ) =F  nc  T  <  /  )  *K 
f  NC  TK ( 2  >  » 1  , /RS 
Fi«iCJR<*0=TZ/KS 
FimCTK(6)«//KS 
F  UC  I  K  (  8  )  «F  NC  TK  (  6  )  *1  Z 
FiiOlT  (31  =  1  ,/RAO/Is 
F  NCl  1  (H)«H/RAD/l!» 

FwCTT(7)=Z/RAD/Tb 
f  NC II  (  8 ) »F  NCTT ( J ) »k 
F  NC  T  Z  ( 5 ) » | , /ZS 
F  N  CUt6)»K/zS 
f  .M  C  1  Z  (  7  )  ■  T  Z  /  Z  S 
FnCTZ(0)»FNCTZ(/)*K 
l»0  1  1  •  I  ,  0 
M (  *  , l  )  *FnC  T  m  (  1  ) 

b  (  2  ,  I  )  *FnC  T (  1  ) 

U  ,  1  ♦  b  )  «FnC  T  T  (  1  ) 
i)  (  3  ,  I  ♦  1  A  )  «FnC  T  Z  (  1  ) 
b  (  S  ,  1  )  *  FN C  T  T (  1  ) 

0  (  *1  ,  1  ♦«  >  bFuCTK  (  I  J-FNCT  (  I  ) 
b(b|i)  - F  N  C  T  Z (  1  ) 

b(!i,IMfc>=FNCTK(lj 
b  ( t ,  1  ♦  8  )  *  F  i*  C  7  /  i  I  ) 

1  b(o,  1*I6)*FNCTT{  1  ) 

IF  C  i  SO  H  3,  1  3 , 6 
13  PR»PD(KK,2) 

Dj=.b«PO ( KK  ,  1  ) / (  1  ,  ♦  pH  ) 

Dm«0j»2./(  i  .-2,*Pk» 

0  I  *  o  H  •  (  1  .  -  p  |<  ) 

0  2  *  0  M  *  P  R 
t>0  2  I  ■  1  i  2M 

bU(I,l)»Dl*(j(  (  ,1  (♦0  2*(b(2,()'»B(  7  »  (  )  ) 
Hb(Z,()*Dl*n(2,l)*U2MH(l,n*t3t3iI!! 
bblJ,I)*Di*n(j,j)>02*UMl,l)>H(2»l)) 


I 
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b6  . 
b  7  • 
bu. 
bV  . 
6U  . 
6  l  . 
62  . 
6  J  . 
64  . 
6  b  » 
6  6  . 
6  7  . 
68  . 

6  V  . 

7(j. 

7  1  . 
72. 
/  J  • 
7  4  . 
7b. 
76. 
77  . 
7  a  • 

7  V  . 
80  . 

8  1  . 
82  . 
8  J. 
84  . 
8b. 
8  6. 
8  7  . 
8b. 
8  V  . 
VO. 
VI  . 

V  2  . 
Vi. 

V  4  , 
Vb  . 

V  6  • 

V  7  , 

V  8  , 
VV  . 

lOU. 
10  1. 
102. 

1  Ui  . 
104. 
10b. 
106. 
107. 

1  0  8. 

1  OV  . 
110. 
111. 
112. 


8  8  1 

4  , 

1 

) 

*03* 

8(4,1 

) 

08 

< 

b  , 

1 

) 

=  0  i  • 

8  l  b  ,  i 

) 

bn 

i 

0  . 

1 

) 

=  0  3  • 

8  (  6  ,  ] 

i 

(i(j 

1  u 

b 

U 

( i  t  it 

.  f-  0  . 

1  )  1)0 

1  i  7 

l\  r 

A 

*  " 

1 

0 

■  •  j  • 

1  z 

I  T 

( 

1 . 

1 

) 

“COG 

(  A2  )  • 

1  M  (  A  ) 

r  r 

( 

2  , 

1 

) 

-  -  C  0 

S  <  A  2  ) 

• 

CUb  (  A  ) 

t  r 

» 

i  . 

1 

) 

3  “  5  1  N  (  A  2  ) 

Ml  I  .2)»C05(A3).Cub(A)^blNlAi)*SlN(A2)*ilN(A) 
Ml2,.*)«t0b(»O)*biN(A)-siN(Ai)*blN(A2)*C0S<A) 

I  I  ii,2)«biN(A3)*CUS(A2) 

TT(|,J)aC0StA3)*SiN(A2)*SlN<A)-SiN<A3)*cC»biA) 

'  1  (2i3)--COb<Ai)*S|M(A2)*CuS<A)-blN<Ai)*bl^(A) 

1  f(i.3)-C0S(Ai)»Cub(A2) 

>  1  <  1  •  I  )  ■  T  I  (  1  ,  1)  •  1  I  (  1  ,  1  ) 

r  <  i ,  2  ) « 1 1  u  ,  2 ) « r  r  t  i ,  2 ) 

l<l,J)*TT<l,J).M(l,J) 
T(1,4»«2,*tf(1,i).it(i,2I 
Jllib)»2.*Tl(l,l)*TT(l,Jl 
I  (  1  i  6  )  «  2  •  •  T  I  1  1  ,2)  *T  I  (  1  ,  J) 

r  t  2  ,  1  )  -  1  f  (  2  ,  1  )  •  1  I  (  2  ,  1  ) 

f(2,2)«ri (2,2)«TT(2,2) 

M2,3)»rr(2,i)*lT(2,3) 

U71'l].2.«Ttl2,l)»lT(2,2) 

H2,!il«^..TH2,l).f1(2iJ) 

T  I  2  , 6  )  »  2  •  •  T  1  <  2  i  2  *  »  1  M  2  ,  J  ) 
J(i,l)«rT(i,i).f],j,i) 

T (  i  i 2 ) *  T 1  (J|2)*lT<3|2) 

T  (  J  ,  J  )  *  T  I  (3,j)*T  T (  j  ,  3  ) 

T  (  3 ,4 )«2. «TT ( 3 , i ) •TT ( 3,2> 
l(J,b)*2»*TT  (  3  ,  1  ) • T  T  (3,3) 

T  (  3 , 6  )  *2.*  T  T  ( J ,2  »  *T  1  (  3 , j) 

T  ( *<  ,  1  )  (  l  ,  1  )  •  T  I  U  ,  1  ) 

M  4 , 2  )»T  T  (  i  , 2  )  •  T  1  (  2 ,2  ) 

1  t  4 , 3  )  «  T  1 1  l ,3)*T> (2,3) 

T(4,4)«Tl(|,2)*TT(2,l)*TT(l,l)*TT(2,2) 

T(4,b)«TT(j,J)*Tl(2,l)'»TT(l,l)*rl(2i3) 

1  (  4  i  <»  )  *  T  T  (  i  »  3  )  •  T  T  (  2  »  2  )  ♦  T  T  (  i  ,  2  )  •  T  T  (2,3) 
Mb,l)»U(i,i)*Tr(j,l) 
f  (  b  ,  2  )  ■  I  I  (  i  ,  2  )  •  T  r  (  3 , 2  ) 

1  (  b  »  3 ) *  T  T (  1  , J)*TT  (i,J) 

T(b,4)»Tr(i,2)*TT(3,l)*TT(l,l)*TT(i,2) 

1  (  b  »  b  )  »  T  1  (  i  ,  J  )  •  T  r ( 3 , l ) ♦ T  3  (  l  ,  1  |  • T  T  (3,3) 

T ( b  ,  6 ) ■ T  T (  i  • J ) • t  T (  3 , 2  )  ♦  T  T ( 1 , 2 ) • T  T (3,3) 
T(6,1)»1T(2,1)*T  1(3,1) 

T  (  6 , 2  )  *  T  T ( 2 • 2 ) • T  T (3,2) 

I  ( fc , J ) «T  r ( 2 , 3  )  *T  1  (  3 , 3  ) 
T(6,4)-Tl(2f2)*Tl(3,l)*iT(2,l)*TT(3,2) 
I(6,b)»Tr(2,3)*ll(3,l)*TT(2,l)*TT(3,3) 
T(6,6)«TT(2,3).TT(3,2)*lT(2,2)*Tl(3,3) 

0UM  =  L2*L3-ti*Li*8l*(Jl-t2*L2«P3*H3 
I”L)*L2*(2.*8l*H2#PJ*P2*p2) 

0  t  1  »1  )*Ll*£.**lLJ-PJ*P3*t2)/UOM 
0  (  1  ,2  )  »L  1  «t2*  (  P2*P  J*t 2* Pi  *t.3  )  /OUM 
0(  1  »  J  )  ■  t  1  *L2*L  3*  (  PI  •P3*P2)/1)Um 
l»(2,2)»L2*C2*(LJ-P2*P2*L1)/0Um 
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1  3. 

0(2,3)=t2*t3*lPl»P2«tl*P3»t2)/UUf1 

1  M  . 

Ulj,j)»t3«£3»U.2*Hi  •  P  1  *  L  1  l/DUM 

1  b  . 

l>(2,I)=0<!,2> 

1  6  . 

U  1  3  ,  1  ) -0 (  l  , 3  » 

1  7  . 

U(J,2)=D12,3> 

1  6  . 

0  1  i  ,  b ) *0 1 

1  b  . 

0  l  b  .  b  )  =  &  2 

20 . 

0  l  t>  i  6  )  =  0  J 

2  1  . 

00  10  11 

22  . 

3 

U(  1  ,  C  )  » 1  . / Kb 

2  3, 

H  (  i 

2b  . 

b  (  2  ,  1  )  =  1  .  /  k  A  i 

2b. 

t!l2,2)=R/NA0 

2  6  . 

li(2,3)*W/KAD 

2  /  . 

8  <  * ,M)*B<2,3>*K 

2  6  . 

8  1  2  ,  /  1  ■  1  •  /  (  k  A  0  •  1  b  ) 

2  V  . 

0i^»6)“n>2|  /  )  •  k 

3  0  . 

IU  3 , 3  )  =  o  (  2  ,  /  1 

3  i  . 

b  l  3  ,  M ) *  0 (2,6) 

3  ^  • 

ii  {  j  ,  b  )  «-b  (  2  ,  1  ) 

33. 

8  l  3 i 6 ) *  1  . /kb-K/h  A  o 

3  *i  • 

1>(j,/)s-B(2,3) 

3b. 

l>(j,t;)«=B(3,6l*T2 

3  6. 

If  t  I 50 ) b , 8  ,  V 

37  . 

8 

PH«F’|3  «KK  ,2  > 

3  fc  . 

UUM*HOUK  ,  1  )  •  h  /  (  1  ,  -  PK*Pk 1 

3  V  . 

80  M  1=1,8 

MU. 

11  b  (  1  ,1  )=UUM*(rt(  1  ,1  )  ♦  P  R  *  t'  1  2  »  1  )  ) 

M  1  . 

Ul3l2,I  )=UoM»(H(2*l  )  ♦  P  K  •  b  l  1  ,1  1) 

M2  • 

M 

8  b  1  3,  1  )  «  U  U  M •  .  b  •  (  1  .  - P H ) • b  <  3  »  1  > 

M3  . 

00  TO  b 

Mb  . 

V 

If  iM<  .  L(j  .  1)00  TO  10 

Mb  . 

A=A|-TU-TS*T2 

Mfc  . 

M  (  1  •  1  )  *  5  1  N  (  A  ) 

/  . 

t  1  l  1  » 2  )  »COb ( A ) 

Mb  . 

M  12,1  1  *”CUb 1  A | 

M  V  » 

,Tl2,2)=bJ^(A) 

b  U  • 

10 

M  1  ,  1  )  =  T  f  (  l  ,  1  )  •  j  I  (i  ,  l  ) 

b  1  . 

1  t  1  ,2  1  *  Tl  i  l  ,  2  )  •  T  T  (  l  ,2  ) 

b  2  » 

I  (  1  ,  3  1  «2«  •  T  T  <  1  ,  1  >  •  1  1  (1  ,2  > 

b  3  « 

U2,l)»Tll2,l)«Tl(2,l) 

bM  . 

1  (2,21  =  1  1  l 2 , 2  )  • T  T  (2,2) 

bb. 

T(2,3)«2»*T1<2,1>»TT<2,2> 

b6  » 

•l3,l)=lt(i,l)»iT(2,l) 

57  . 

M3,2)  =  Tf(i,2)*Tl(2,2) 

bb. 

T(3,3)«TT(i,2)«tT(2,|)mT(1,1)»TI 

b? . 

0  U  >'i  =  L  2  -  P  1  ♦  P  1  •  t  1 

60. 

0(  1  ,1  ) ■£ 1 • £  2  2  U  U  M 

61  . 

0 (  1  ,2 )«P  )  *o  i  1  ,  1  ) 

6  2. 

0  *  2  ,  1  ) -0  »  1  ,  2  > 

6  3. 

U12,2)«E2*E2/DUP. 

6  M  . 

0  (  3 , 3  )  »  0  1 

6  b  * 

1 1 

ou  lb  1 ■  1  ,  1  b  1  N 

6  6. 

0  (.  lb  J  a  1  ,  'j  A 

67. 

i>p  i :  ,j » -o. 

6  U  « 

Uu  1 7  K» 1  ,  1 D  1  M 

6  V  . 

1  / 

UP(  1  ,  J ) »  0  P (  1  ,  J)*U(  1  ,K)*1 ( j,k  ) 

*2) 


1  7  0  . 

i  e 

com  i nul 

1  /  1  . 

1  0  *  1  0  1  M  ♦  1 

17-!. 

00  20  1  ■  1  0  ,  r.  X 

1  7j. 

01/  20  J  ■  1  i  N  a 

174. 

2u 

op  u  ,  j  )  *0  (  1  ,  1  )  •  T  (  J  t  j  f 

1  7b. 

00  22  1=1 tNA 

17  6. 

1'0  22  J«  1  ,  nx 

177. 

0  (  1  ,  J  )  •  0  * 

1  7  b. 

Ou  21  1 , NX 

179. 

2  1 

0<  1  iJ)«Ut  1  ,  J  )  ♦  T  (  1  ,  K  )  •  D  P  (  N  |  J  ) 

18  0. 

CONT INUfc 

1  fa  1  . 

00  2  3  1  *  1  , Nx 1 

18  2. 

1  1  ■  1  ♦  1 

183. 

Oo  23  J« 1 i ,nX 

184. 

23 

0( J,  l )«0(  1  , j) 

1  bb. 

00  2b  1  *  1  , NX 

1  86  . 

Oo  2b  Ja 1  , n  Y 

187. 

bH(  1  ,  J)«0. 

18b. 

Oo  24  K*1 ,NX 

189. 

24 

Obt  1  ,  Jlabbl  l  1  ,K)*b(K,J) 

1  VO  . 

2b 

Con i i nul 

191. 

b 

8  L  1  0  K  N 

1  V2. 

t  NO 

fc  NO  Of 

compilation:  no  diagnostics. 
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!  . 
7* 
3 . 

*4  • 

5  . 

6  • 
7  . 
A  . 
9  . 

I  n. 
t  t . 
i  ?  . 

1 3 . 


c 

c 

c 

c 


SUBROUT  I  Hi  rJTrG!  |  OIM.KK.  NO!  ,  NDF.SU 
for<M5  stiffness  Matrix  of  elements 


0,'MFf  SION  51.  (21»?1),RP(6,j><|)|AA(2),t,'X(21|?1) 
COMPPN/TANAl  /  RAfl.RS  ,TO  .TS.ZS.NX  ,HNB 
coMHON/puTH/Lor (*snon> .nTnn.nb.np.minuh 

C  nMMi.il/RG  AMD/(.  (fl.Ml  ,R(6.2l)  ,H 

0* ta  AA/. 5773502691 fl 96 2 6  ,  - . 5 7 7 350269 l 8 96 26 / 

CAl L  C0ORD3INDI ,R| , Tj , Z 1 ) 

CALL  C0ORD3 ( NOT , R? , T2 , 72 ) 

RS*.S* (R7-R1 | 


I  1  . 

R0» , 5* ! R2*R  1  1 

1  S  » 

IF  ( T  2  . L  T .  T|  IT?*  6. 28318330' 

1  A  . 

T  S  =  «  5 • ( T?-T  |  1 

1  7  . 

T  0  « . 5 • ( T?4T  1  I 

1  «  . 

7 S * . 5*177-7)  j 

1  9 . 

If  1  1 0  I M  .ro.  7  I 7Sa  )  . 

20  . 

M  C  X  ■  1  l)  I  M  -  | 

?l  . 

NP I  ■NP  -  | 

?7. 

0  0  7  I  a  1  , M P 

2  3; 

0  0  2  J  » I  , N P 

2*. 

2 

SL 1  1  .  JIbO, 

25  . 

On  7  I«) ,2 

?6. 

R  =  A  A  (  1  1 

27  . 

RAPpRU*PR*r 

2fl. 

nunieRAn*Rs*Tri*zs 

2  9. 

no  6  J* 1  , 2 

30 . 

T  *  A  A  (  J  ) 

1 1  . 

00  5  K ■ I  , N C  X 

37. 

Z  ■  A  A  (  K  ) 

33. 

f  Al  L  MTRXP(  IDIM.Klf  ,R,TfZfRR| 

3  1  , 

00  1  NRs 1 , MP 

35. 

00  1  NCrNR.NH 

36. 

D  U  M  ?  ■  0  , 

37. 

0  0  3  N  a  |  ,  n  y. 

3a. 

3 

OllM7*DUM?*r.  (N.NRI  •rntN.NC  1 

39  . 

1 

SL  (  NR.Nt  )*SL  (MP  ,nc  )  ♦PL'MI  *DUM2 

in , 

S 

CONTlNUr 

1 1  . 

6 

comt  r  nuf 

12. 

7 

CONTI  HUE 

*M. 

00  «  I  B I  , NP I 

11 . 

J*  r  «•  i 

15. 

oo  e  k»j,n p 

16. 

H 

SL(K  .  I  )«SL( I  ,K  » 

17  . 

I  1  ■  -  1 

in . 

oo  i  i  i « i  , up 

19  . 

1 1 » 1 1  ♦  i 

So , 

1 ?• (  1  • I  )  / | 0  I M« | 

5  1  . 

00  10  J  B I  , N  P 

S  7  . 

OX  (  I  ,  J  I  .0. 

S3. 

00  9  6B I , MR 

Si . 

< 1 3< ♦  1  1 *NR 

S5  • 

0X(  I  ,  J) bDX ( I ,J)*C(K,I2)*Sl(K!  , 

J) 


1 


r 


5A, 

OKI  I  *  J 1 s  0  x  (  I  , JI/HMH 

S  7  . 

1  n 

C  n  N  T  I  N  u  f. 

SO. 

IF  (II  .  F  Q  ,  N C  X  I  1  1  -  -  1 

S<». 

1  1 

CONTI NUr 

AO. 

1  1  «=-  1 

6 1 . 

0()  M  1  »  1  ,  >!P 

A?. 

1  1  «=  1  1  ♦  1 

A7. 

!?*(  l-i  )/min*i 

AH  . 

0  0  17  J  3 1  , N  P 

AS. 

St  1  J  , l  1  »n. 

A  A  * 

00  | 2  Ka 1 

A  7  . 

K  1  =  «♦  1  1  •  N  fl 

An  . 

I  7 

SHJil  l«SLIJ,|  l+DXIJ, 

60  . 

Sl  <  J  .  1  1  aSL ( J .  1  ) /HNR 

7  n , 

1  7 

CONTI NUF 

7  1  . 

IF  (II  . F  Q .  N  f X )  1  1  a  - | 

7?. 

1  M 

CO NT  I NUf 

73  . 

00  IS  1  r  1  ,  HP  l 

7^. 

J  ■  1  ♦  1 

7  G  , 

0  0  IS  K  r  J , N p 

7  A  . 

1  s 

Sl  <  1  .  K  I  aSL  (  tf  1  l  1 

77. 

W  F  T  1 )  W  N 

7  n . 

FN  I) 

fno  nr  compilation;  mo  diagnostics. 
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1 . 

SUBROUTINE  STwtsSlKK.NOD.STRS) 

2. 

COMMON/ NOth/LOC (2&000) , NTNN , NB , N* , M | Num 

3. 

COMMON/T*NAl/HOt«S,TO,iS.ZS,Nr,HNB 

S. 

COMMOn/JMFER/NH.Ns.nZ.NwOP, | dim 

b. 

common/ AO*Ln/D( 3.600) 

6. 

COMMOn/FUnk Y/NBkC (500) ,nRE0(500) . NOE* i bOO) 

7. 

common/ 8<»  and /cu*8),o<6, 21), h 

8. 

DIMENSION  N00(8)  *DSPI21)  ,A(21),BB|8,2i),STRS(*) 

9. 

C 

EVALUATES  COMPONENTS  of  STRESS  along  CYLINDRICAL 

10. 

1  1-0 

I  1  . 

Do  11  J- 1 'NO 

12. 

NDO-NOO ( j ) 

13. 

MTT=NTNN"N0ER(NDD)41 

It, 

DO  10  K-l  , I D 1 M 

lb. 

I  1  ■  1  Ul 

16. 

10 

DSP (  U  )  -D ( « ,M  T  T } 

17, 

1  l 

CONT InUE 

10. 

I  1  --  1 

19. 

DO  IS  J-l.NX 

20. 

A  (  J  1  -0. 

21  . 

Jj-< J-l  )/lOlM*l 

22. 

I  I  -  I  I  ♦  A 

23. 

00  13  K-  l  ,NB 

21. 

13 

A(J)-A(J>*c<JJ,K>*DSP(ll*NB*K) 

2  S. 

IF  (11  . Lq ,  ( IDlM-l ) ) I l.-j 

26. 

lb 

CONTINUE 

27. 

call  mtr*b< ioim.k* ,o. ,0. ,0. ,bb> 

20. 

DO  19  J-l.NY 

29. 

STRS( Jl-O. 

30. 

DO  18  K-l  , NX 

31  . 

10 

STRS(J)-STRS(J)4BB(jtK)»A(K) 

32. 

19 

STRS( J1-STRS1 J)/HnB 

33. 

return 

31. 

end 

END  or  COMPILATION: 


NO  DIAGNOSTICS. 
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1 . 

2. 

c 

SUBROUTINE  PSTNESf 101M,hMa,B,P1 

3. 

c 

computes  magnitudes  and  direction  c 

9. 

c 

of  principal  stresses 

5. 

c 

6. 

CUMM0n/nAnaY/CH(3,3),NA 

7  . 

DIMENSION  8(6) , p ( 3 ) 

a. 

If  <  ID1M  .EQ.  2 ) G o  TO  13 

9. 

A  1  »B ( 1 ) *6 ( 2 ) *B  (  3  1 

10. 

1 1 . 

lB(9)*B(9)-B(5)*B<5)«B(6)*B(6l 

1  2 . 

A3»B<i)*fl<2)*B(3U2.*BUl*BC»)*B<6>. 

I  3, 

1“B(2)*B(S)#b(S)„B(3)*B(s)*B(9) 

19. 

P  (  3  )  ■  a  1  /  3  , 

IS. 

OM«P ( 3 1 

16. 

00  S  1 ■ 1  •  9 

1  7. 

FX*P(3)*IA2*P(3)*(P(3)-a11)-A3 

IB. 

FPX-A2*P<3)*(3.*P(3)-2.*AI> 

19. 

P(3)-P<3»-FX/FPx 

20. 

5  CONTINUE 

21  . 

PRINT  6 , U  m , p ( J )  ,px 

22. 

6  FORMAT  ('  FIRST  P(3)  «  *,E1Q.5,*  •• 

23. 

1 E 1 0 . 5  ,  •  ••  ITER,  ERROR  ■  ' • E 1 0 . 5 ) 

29. 

IF  ( AbS ( F  X )  . L 1 »  1 50 *  ) GO  TO  8 

2S. 

Om-P ( 3 ) 

26. 

00  7  1  ■)  ,9 

27. 

FX-P(31*(A2*P(3)*(P|31-A1))-A3 

2B. 

FPX-A2*P<3).(3..P(3!-2,*A1  ) 

29. 

P ( 3 ) "p ( 3  1  -FX/F PX 

30. 

7  CONTINUE 

31  , 

PRINT  6  ,  Dm , P ( 3 1  , F  x 

32. 

8  B lap ( 3). A  I 

33. 

B 2 ■ P ( 3  )  1  *A2 

39. 

OMaSQRT (B|»b1-9,*B2) 

35. 

P  <  1  » ■ . 5*  < -8  l  -DM ) 

36. 

P ( 2 1  * . 5« ( -B  1  *DM ) 

37. 

IF  (P(3)  , G  T .  P ( 2 )  ) GO  TO  10 

33. 

B  2  ■  P  (  2  > 

39. 

P (2)«P(3> 

90, 

P ( 3 1 "82 

9  1  . 

IF  (P(2)  ,GT.  P< l ) )G0  TO  10 

92. 

B  1  ■  P  <  l  ) 

93. 

P ( 1 1 «P ( 2 I 

99. 

P ( 2 1 "8  1 

95  , 

10  IF  (Mma  • NE «  1)G0  TO  20 

96. 

c 

CR  •  direction  cosines 

97. 

00  12  la!, 3 

98, 

CC ■ ( P ( I )"B(2l 1 *(P ( 1 )-b( 3) ) >B ( 6 1 *B ( 6 ) 

99, 

BBa (P (  1  1-811 > 1 *(P( 1)-B ( 2  »  >-B ( 9) *8(9) 

SO, 

IF  (AbS(bB)  .gt,  abs ( CC ) >  go  to  li 

bit 

AA"(P(  1  )  '8(31  )*B(9)*B(b)*B(6) 

52. 

Bb«(P( 1 )-B(2) ) *B ( 5 ) ♦B ( 9 ) *B ( 6 ) 

S3 , 

GO  TO  V 

59 , 

11  CC*(P(  1  1-8(2) 1 *8(5) *8(91*8(6) 

b  b  • 

A  A* ( P (  1  )“8( 1 1  )*B(6l*B(9)*B(5) 

OSINES 


final  p(3, 


S6. 

9 

57. 

5B  , 

59. 

60. 

12 

61  . 

62. 

1  3 

63. 

68. 

65. 

6  6. 

67. 

68. 

69. 

70. 

71  . 

72. 

73. 

78. 

lb 

7b. 

76. 

77. 

I  8 

78. 

20 

79. 

0m«5QrT(Cc*CC*B8*BB*AA*aA) 

CR  (  1  ,  |  )  »Cc/DM 
CR ( 2  %  I  )  «Aa/OM 
CR ( 3  i  |  ) -Bb/DM 
Con  t i nue 
0  0  TO  2  o 

dm. SORT ( • 25»  <  B  < | >-6(2) )»*2+B(3)«B(3)) 
1  >*.S*<B( I  )  ♦B ( 2 )  )*DM 
P(2>".S»*BC l  )  *8 ( 2 ) WDM 
1 F  ( MMA  • NE  «  1 ) GO  TO  2Q 
DO  1*4  1.1,2 
A  A  - P  (  |  ) -B (  i  ) 

BB«P ( I l -8 ( 2  ) 

IF  (ABS(BB)  ,GT.  ABS(AAj)G0  TO  15 
DM.SQRT (Aa*AA*B(3)«b( 3)  ) 

CR ( 2 , 1  )  «Aa/0M 

CR(  1  » I  1  «B ( 3  )  /dm 
GO  TO  H 

Dm-SQrT (Bb«BB4B(3)«B(3)  ) 

CK (  1  ,  l  1  .bb/dm 
C  R ( 2  ,  |  ).6(3)/dm 
C  ON  T I  NO> 

RETURN 

End 


end  of  compilation; 


NO  DIAGNOSTICS, 
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I  . 
2. 
3. 

<4  • 

5. 
A  • 
7. 

n . 

o. 
10. 
1 1 . 
12. 
13. 
1  A . 
IS. 
1  6. 
17. 

1  A. 

19. 

20. 
21  . 
22. 
23. 

2  A . 
25. 

2  A  . 
27. 
20. 

29. 

30. 
31  . 

32. 

33. 

31. 

35. 

36. 

37. 

3  A  * 
39. 
*♦0. 
A  1  . 
A2. 
A3. 
AA. 
A  5  . 
A  6  « 
A  7  . 
A8. 
A  9  , 
50. 
51  . 

52. 

53. 
5  A . 
55. 


SUBROUTINE  FACTOR 

C  •••••••••••••••••••• 

C  MATRIX  FACTORIZATION 


COMMON/fUNKY/NR0W(500)  ,NReO(500>  ,N0eR(500) 

common/jnfer/nr,ns,nzx,nwop,  I  DIM 

C0MM0N/A8ALN/DI 3,600) 

common/ janet/nspd.nwsdu 00)  ,SPO(  100,3  1 

COMmON/BOPaP/NLOD' 100) *  C  NL  0 ( 1001 , NE  ,NER 
C0MM0N/RUTH/U0C( 25000J , N TNN . NB  ,  NP  ,  M 1  ST Y 
DIMENSION  A ( 3  i  300 )  .1(3,3) ,C<3,3)  ,JlT(A) ,LC0N(8)  , 

I  NODCfl) ,NSPT (500) ,S(3,300) , SL ( 2 A . 2 A ) 

EQUIVALENCE  (DM, 301), S) 

REWIND  15 

rewind  ia 

RFWJND  13 
REWIND  12 
REWIND  11 
REWIND  10 
J  I  T  (  1  )  »  0 
J  I  T ( 2  )«0 
J | T ( 3 ) «0 

M  1  M  |  «0 
•CD-0 
NE  E  •  1 
MEE«1 
N  Z  *  0 

DO  70  1 «  1  ,  NTNN 
N  N  •  0 

IF  (I  , F  Q ,  NTNNIGC  TO  6 
NN.-NROW ( J ) 

NROW ( 1  J  «NZ+ 1 
NSPT (  1  I «NZ+  1 
NZ«NZ*Nn> 1 

6  NN3« ( NN*  1  )  • J  0 1 M* 1 
DO  7  J.  J  ,  1 0  1  M 

DO  7  K. I  ,NN3 

7  S ( J , K ) «0  » 

JaNREO (  I  1 
PRINT  100,1 , J 

100  FORMAT  (6H0N00E  •  I  A  ,  J  1  H  (ORIG*  NO,,fA,lH)) 

CALL  ELTONO  (  1  D I M , j , LCON ) 

PRINT  I  01  , (LCON(  I  1  )  ,  If. J  ,NR  ) 

101  FORMAT  (3AH  ELEMENTS  CONNECTED  TO  NODE  --  ,8(lA,lH,)) 
DO  1 5  K.  1  ,  NB 

IF  (LCON(K)  *EQ«  0)50  TO  15 
KK-LCONIK ) 

CALL  NODES  ( IDIM.KK  ,NOD) 

PRINT  102.KK, (NOD( I  1 )  , l  I-|  .MB) 

102  NDUNoo!fIH  N00F:S  of  element  ,ia,ai.  ,«ma,1h,)) 

NOF.NOD ( NB- 1 ) 

I F  ( NER  , EQ .  1 ) GO  TO  56 

K S • ( K K - j  ) /NER*  I 

IF  (KS  , EQ ,  K  D  )  GO  TO  57 
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56. 

NUMb(KS-1  ) /M 1  ST  Y ♦  1 

57. 

1DAUM«NUM-»12 

5  A  . 

jNT«K‘J-<NUM-l)*MlSTY-JITtNUM) 

5". 

399 

FORMAT  (•01DRUM»»i13,»**»JNT«*.13) 

60, 

PRINT  399 , f DRUM , JNT 

^  1  • 

JlT<NUM!r-KS-<NUM-l  )  *MI ST Y 

67. 

IF  (JNT  .LE.  0 ) GO  TO  Si 

63  . 

IF  (JNT  .  E  0  »  1  ) G  0  TO  5  M 

6m  . 

JNT  «  JNT -  1 

65. 

DO  50  1  A ■ 1  ,  JN  T 

6  6. 

READ  ( 1 0  RUM  1  (  (SL(L1 . L  2 ) » L 1 » l , NP ) , L  7*  1  . 

67  . 

50 

CONTI NUF 

6  A  • 

GO  TO  SM 

69  . 

51 

JMT«  1  - JNT 

7n , 

DO  S  2  1  A  *  1  , J  M  T 

71. 

BACKSPACE  IDRUM 

7  2, 

57 

CONTI  NUF 

73. 

5M 

READ  (  IDRUM)  <  (SL(L1  . L2  )  ,L 1 • 1  i  N  P )  ,L2"1  i 

7  A  , 

((D«kS 

7  5, 

GO  TO  57 

76. 

56 

Call  INTEGI  IDIMjKK.NOI  ,ndf,sl» 

77  . 

57 

DO  A  M«|  ,NR 

7  A  « 

IF  ( J  «EQ.  NOD(M)  1  GO  TO  9 

79. 

A 

continue 

An . 

9 

ND I  A  »M 

A  1  . 

DO  10  KM« 1 , 1 D 1 M 

82. 

K  N ■ ( ND 1  A- 1  ) • 1 D I  M  ♦  K  M 

83, 

00  in  L  M • 1  ,  1D1M 

AM  . 

L  N ■ ( N  D 1  A “ 1  )*1DIM*LM 

AS  , 

in 

S«KM,LM).S(KM,LM14SL1KN,LN) 

6  6. 

IF  INN  ,EQ,  0 1  GO  TO  IS 

87, 

DO  18  M. 1 ,NB 

8  A  . 

NODM»NOD ( M ) 

89, 

IF  (1  ,GE.  NOER INODM)  )G0  to  im 

90, 

NRWrNROA (  1  )  -  I 

91  . 

DO  13  L« 1  ,NN 

92. 

NRW-NRW* 1 

93. 

IF  (  L  0  C ( N  R  W )  ,NE.  NOER I NODM )  ) GO  To  13 

9  M  « 

DO  12  K  M ■ I  ,  1 D I M 

9S  . 

KN« ( ND 1  A- 1  ) • I 0 1 M*KM 

96. 

DO  12  MM» 1 , I D 1 M 

97  . 

LN«(M-1  )  • J  D I M4MM 

9  A  , 

L  M  -L  • ID  1 M ♦ M  M 

99, 

I  2 

S(KM,LM)»S(KM,LM)*SL(KN,LN) 

100, 

GO  TO  1M 

10  1. 

l  3 

CONT INUE 

102. 

1  M 

continue 

103. 

15 

CONTI NUC 

1  0M  , 

IF  ( J  »NE.  NLODINEE  )  1G0  To  16 

i  ns. 

S (  1  ,NN3 J -CNLD (NEE  ) 

106. 

NEE-NEE* 1 

107. 

16 

WRITE  (  1 2 ) NN  3 • (  ( S  <  L  1  ,  L  2  )  iLI"l  •  I D  I  M  )  ,  L  2 

1  0  A  » 

IF  (MEE  .GT.  NSPD1G0  TO  38 

I  09. 

IF  (J  .NE.  NWSD(MEF) )G0  To  18 

110, 

DO  17  K  a  1 , I 0 1 M 

111. 

IF  (SPD(MEF.<)  *GE.  100. )G0  TO  17 

112. 

DO  26  L- 1 , 101 M 

113. 

S(L,NN3>»S(L.NN3)«.S(L,K|*SPD(MEE,K) 

1  1  9  . 

S  (  K  ,  U  ■  0  . 

I  I  5  • 

26 

5 ( L  *  K 1 »n  • 

IU, 

S  (  K  ,K  »  ■  1  . 

I  I  7. 

S  (  K  ,NN3 1 -SPD ( MEE ,K 1 

tin. 

I  7 

CONTI NUE 

119, 

MEE.HEE* 1 

1  20. 

IF  (MEE  .GT.  NSPOlGO  TO  30 

12  1. 

I  8 

MU-MEE 

I  22. 

NU*  1 

123. 

NRW  1  bNKO'.V  (  1  )-l 

129. 

00  28  K-MU.NSPD 

1  2  S  • 

NWS-NWSf)  (  X  ) 

12A. 

DO  30  LoNU.NN 

127. 

NRWbNRW|*|_ 

1  28  . 

IF  ( LOC ( NR w )  .GT.  NOER  1  NWS ) 1  GO  TO  2 8 

I  29. 

IF  (LOC(NRW)  .'IE,  NOER1NWS11GO  TO  30 

I  30. 

MM-U 

13  1. 

M3«(NRW-NR0W( 11*1) • 1 0 1 M 

132. 

00  39  I  A  -  I  ,  1 D  I  M 

133. 

IF  ( SPD ( K  .  1  A  1  .GE.  100.160  TO  39 

139. 

00  36  1  R.1  , I 0  I  M 

1  35. 

S(!P,NN3)«S(IB,NN31-5(IB.M3*!A)*SP0(K,U) 

13a. 

IF  (S(TB,  101-1.  .GT.  .001  IS!  1B,M3*JA)*0. 

137. 

36 

CONT I NJE 

138. 

39 

CONT I NUf 

139. 

GO  TO  37 

I  90. 

30 

continue 

I9J  . 

GO  TO  38 

192, 

37 

NU-MM+ 1 

193. 

28 

CONTI NUE 

199, 

3  8 

IF  (  M  1  M i  . EQ .  0 1  GO  TO  62 

1  95, 

JILL-1 

196. 

1 1  - :  - 1 

I  97. 

DO  60  K> 1 , I I 

I  9  fl  , 

RE  AO  (  1 0  I  M  N  1  .  (  (A(L1  »L2  1  »L  1  •  1  .  I D I M 1  ,  L2»  1  |NN|  1 

1  99. 

1  1  (Bill » L  2  1  ,  L  1  •  1  •  I 0 I M  )  ,L2*l , I  0  I M 1 ,0ET 

160, 

ispt-nspt  <K  > 

16  1, 

IF  (LOC (  1  SPT  1  .NE.  I  1G0  TO  60 

162. 

lTFR-CNSPT(K)-NRPwmcn-IDlM 

163, 

NS  P  T ( K  1  bNSFT ( K  1  ♦ 1 

159. 

IF  (JILL  .GT,  NSPDJGO  TO  92 

155, 

DO  39  LI-JILL.NSPO 

I  6  A  , 

NWO-NWSp ( L  1  1 

157. 

KEE-L  1 

158. 

IF  (NOER(NWD)  .EQ,  K)GO  TO  35 

I  59, 

IF  (NOER(NWD)  »GT,  K)GO  TO  92 

160, 

39 

CONTI NUE 

16  1. 

35 

J  !  LL  *KEE ♦  1 

162, 

DO  9 i  I  A  «  I  , 1 D  1  M 

1  63. 

IF  ( SPO ( KEE i 1  A  1  .GE.  JDO.IGO  TO  91 

169, 

00  90  10- 1  ,  1 D  I  M 

1  65. 

S  (  IB  ,NN3 ) -S ( IB ,NN3 1 -SPD(KEE , 1  A  1 *A ( 1  A , I  TER* IB  1 

1  6A  . 

90 

a  <  i  a  ,  i  ter* 18 1 »o. 

167. 

9  I 

CONTI NUF 

168. 

c 

GAUSSIAN  ELIMINATION 

1  69. 

9? 

DO  21  1  A ■ I  , 1 D  I  M 
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1  70  . 
17  1. 
17?. 
173. 
1  74  . 
I  75, 
I  '  6  • 
177. 
1  7  A  . 
1  79. 
1  «0. 
1  8  1  . 
l  e?. 

183. 
1  04  . 
135. 
1  86. 

187. 

188. 

189. 

190. 
I’l  . 

192. 

193. 
199. 

195. 
1  96. 
197. 

1  9fl, 

199. 

200. 
201  . 
202. 
203. 
209  . 

205. 

206. 
207  . 
208, 
209. 

2  10. 
2  11. 
212. 
2  13. 

219. 

215. 

216. 
217. 
2  18. 
2  19. 

220. 
221. 
222. 
223. 
229  . 

225. 

226. 


I  D *  T TER* I  A 

00  21  1 8 ■ 1  i  ! 0  1  H 

C  (  1  A  ,  l  B  » -0. 

00  19  JOl.lDIM 

19  C(!A,!B)sC(!A.1B>*AUC,!D)*8(!C,IB) 

2!  CONTINUF 

NN0«NN 1  - l TER 
00  29  I  A ■ 1  ,  1 0  1  M 
00  23  18-1, NNO 
0 ( 1  A ,  l  B  1  »0. 

10-1  TER*  IB 

00  22  1C-I.I0IM 

2?  0(|A,!8)-0(IA,IB)*C(IA,1C)*A(IC,I0) 

0(!A|lB)«D(  1  A i I B ) /DET 
23  CONTINUE 
29  CONTINUE 

L  YN-NNO/ I D I M -  l 
IF  (LYN  , EQ •  0)60  TO  32 
NRW-NROW (  I  ) -  1 
J  S  -  1 

DO  31  L« 1 , L YN 
MSPT-  1  SpT  ♦(. 

L  3  «L • IDJM 
DO  2'  M.JS.NN 
NRW«NR6*1 

IF  (LOC(MSPT)  .NE.  L0C(NRW))6O  TO  27 
M  M  ■  M 

M3»(NRW-.NR0W<  1  )  ♦  l  1  •  I  D  I  M 
00  25  I  A  *  I  ,  I D 1 M 
00  25  I R. 1  ,  1 D  1  M 
I L "L  3 ♦  1  R 
1 M»M3* 1 B 

25  S (  1  A  ,  1  M) «S( I  A , I M ) -D< l A  ,  IL ) 

GO  TO  29 
27  C  0  N  T 1 N  U  F 
GO  TO  32 
29  J  S ■ MM ♦  | 

31  CONTINUE 

32  00  33  I  A ■ 1  i  I  0  I  M 
S(!A,NN3'»5(  1A,NN3)*»D(1AiNN0) 

00  33  l 8 ■ I  , I  0  I  M 

1 0“ I  TER ♦ I  8 

A(1A,ID)»C(!8,IA1 /OET 

33  S ( I  A , IB ) -S(  I  A , IB )-0( I  A  ,  IB) 

BACKSPACE  10 

WRITE  (I0)NN1,(<A(L1»L2),L1»1»I0IM),L2«I,NNI), 
1  I  (B(Ll  »L2)  , L I ■  1  i I01M)  ,L2*I . I D I M )  ,  QE  T 
IF  (K  »EQ •  I  I ) GO  TO  60 
BACKSPACE  10 

READ  (  l0)NNl  .  (  (  AILI  iL2)  •(_!*■!  •  I  D  I  M  )  ,L2"I  .NNI  )  , 
I  I  (B(L1  « L  2  )  , L I ■  1  •  IDIM)  ,L2*1  .  I  0 1 H )  »OET 
60  CONTINUE 

6?  OET*S(  1  ,1  )  *S  (  2  i  2  ) -S ( J l2)*S(I|2) 

IF  (  101M  .FO*  3)60  TO  69 
B( 1 . 1  ) - S  C  2 , 2  I 
B ( 2 ,2 )*S  C  1  ,  1  1 
B( I  , 2 ) »  >  S ( 1,2) 
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227  , 

220  i 

2  2  9, 
230  , 
23  1  . 
2  3?. 
2  3  3  . 
23m  . 
2  3  5  , 

236  . 

237  , 
235  , 
239  . 
250  . 
2  M  1  . 
2M?. 
283. 
2  M  M  , 
2  M  5  , 
286  , 
207, 
28R  . 
2  8  9, 

250. 
25J  . 
2  52. 
253. 
2  5m, 
255, 
25*. 
257  . 
2  5  fi  , 

259. 

260. 
26  |  , 
262  , 
263  . 
2  6  M  . 
265. 
266  . 
267  . 
2  6  R  , 

269  , 

270  . 
27  |  , 

272  . 

273  , 
2  7  M  , 

275. 

276. 
277  . 
2  7  R  , 

279. 

280. 
2  R  I  . 
282. 
283  . 


GO  TO  65 

6  «  fl(  I  ,  l  )  ■  5  (  2 , 2  )  •  S  (  3,3)-Sf2*3)*S(2,3) 
fi(1.2>-5M,3>*S(2,3)-Sn,2)*S<3,3) 

8( 1  . 3 ) • S  f  I  ,2)*S(2,3)»S(  I  .  3  )  *5  <  2 , 2 ) 

8  (  2  2  >  -5  (  1  ,  1  )  *S  (  3  ,  3  )-S  (I  ,3)  *SM  ,3  ) 

81?  ,  3  )  *»5  (  !  ,2  )  •  S  C  I  , 3 ) -S (  I  .1  ) *S ( 2 , 3  ) 

B I 3 , 3 ) »0E  T 

0ET«5( 1 , 3 ) • 8 ( l ,3).S(2,3)*B<?,3)*5(3.3)*B(3,3) 

8  |  3 , 2  )  ■  p  (  2 , 3  ) 

R  (  3  ,  1  )  »p  (  I  ,  3  ) 

65  8  (  2, 1  ) -n <  1  , 2  ) 

IF  (I  «EQ.  NTNN ) GO  TO  70 
NN I «  N  N  • 1 0 1 M ♦  1 
N  3  ■  I  0  1  M  *  1 

WRITE  (  1  I  ) N N 1  , <  <S<L! »l2)  ,l  l  *  1  iIDIM)  ,l2»N3.NN3) 
WRITE  M0)NN1,I(S(LJ,L2),L1«1,|DIH),L2»N3,NN3», 
I  (  <  B ( L  I  , L  2 )  , L I  *  1  • 1 0  I M )  ,  L  2 ■  I  ,OET 

IF  <  M | M |  , E  Q ,  1)60  TO  6? 

00  66  K>  1  ,  I 
I SPT  =  NSPT  (<  ) 

IF  (LOC(ISPT)  ,  £  Q  «  1*1)60  TO  68 

66  CONTINUE 
GO  TO  7n 

6  R  M  i  H  1  »  l 

69  REWIND  10 

70  CONTINUE 

c  . . . . 

c  BACK-SUBSTITUTION 

C  . 

00  76  1.1  ,  10  IM 

0  f  I  ,  1  )  -  0  . 

00  7  1  J» 1  ,  1 0  1  M 

71  D(  I  ,  1  )»D(  I  i  1  )*B(  1  ,J)*S(J»NN3) 

76  0 (  I  , 1  ) »D (  1  ,  1  ) / OE  T 
00  75  I  ■  2  ,  N  T  N  N 
NOOE-NTNN-  I  ♦  l 

backspace  10 

READ  I  1 0  >  N  N  1  .  (  t  A  t  L I »L2)  »L I  *  1 iIDIM)  ,L2»I  »NN  I  )  • 

I  < (B(LI  i L  2 )  ,L  1  ■ l i I  0  I M  )  | L  2»  l  ,  1  0  I  M  >  ,DET 
BACKSPACE  10 
00  77  1  A  *  I  , I D  1  M 

0  (  1  A  ,  I  )  ■  0  » 

00  72  1 B  *  I  , I  0  I M 

72  0(lA,I)aD(IA|l)*R(tA,IB)*A(IBtNNl) 

77  0 {  I  A , I  ) .0 (  I  A , I  ) /OE  T 
NRW 1 »NR0W (NODE ) -  1 

nn-nni/jdim 

00  7M  J«  1  ,NN 

j j- ( J - 1 ) •  i  n  l  m 

NRW-NRW  1  ♦  J 
MRW»NTNN-LOC  ( NRW ) ♦ 1 
00  73  I  A ■ 1  , I  0  I  H 
00  73  iRaJ.lOIM 

73  0  (  I  A  i  I  I.OIlAil  )-A<  I  A , J J*  1  B  )  *0 (  I B .MRW ) 

7  M  CONTINUE 
75  CONTINUE 

PRINT  300 


284. 

300 

FORMAT  (*ONODAL  displacements* ) 

285. 

DO  85  J» 1  ,  I  o  I  M 

288, 

PRINT  350, <0< J.K) ,K«NTNN,  1  ,*!> 

287, 

85 

CONTINUE 

288. 

350 

FORMAT  (12EI0.5/U0X, I1E10. S»> 

289, 

RETURN 

290. 

END 

END  or  COMPILATION} 


NO  DIAGNOSTICS. 
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1  . 

2  . 
3. 
B  • 
5. 
A  • 
7. 
B. 
9 . 

10. 
1 1  • 
12. 
13. 
IB. 
IS. 
1  A  • 
17. 
IB. 
19. 

2n. 
21  . 
22. 
23. 
2b  . 

25. 

26. 
27. 
2fl. 

29. 

30. 
3  1  . 

32. 

33, 
3  B  , 
3S. 
3  6. 
37. 
38  . 
39. 
BO. 
B  1  . 
B  2 , 
B  3  , 
B  B  .. 
B5. 
B  6  , 
B  7  , 
Bfl. 
B  9  » 
Bn, 
5  1  * 

52. 

53. 
SB, 
55, 


SUBROUTINE  FC  TOR 

COMMON/FUNKY /NROW (500)  ,NRE0(500) ,N0ER(500) 

COMMON/ JNFER/ NR , NS, NZ ,NWOP,  I  0  I M 

common/ruth/loc i 25000) »ntnn»nb»np, misty 
COMMON/ j a  NET /NSPO.NWSOI  100) ,SPD( 100,31 

COMMON/ ABALN/O ( 3 , 600 ) 

DIMENSION  A ( 3 , 300 ) ,6(3,3)  ,C(3,3) , NSPT ( 500 ) ,$(3,300) 
EQUIVALENCE  10(1,301)  ,5) 

REWIND  10 
REWIND  II 
REWIND  12 
MEE-  1 

PRINT  100, MISTY 

100  FORMAT  MOSTARTINq  point  of  GAUSSIAN  ELIMINATION  IS  R 0 W •  , 
IF  (MJSTY  ,EQ.  1100  TO  9 
MRX»M 1  STY- I 
DO  7  1  ■  |  ,MRX 

IF  ( NR  F  0 1 1  )  .EQ.  NWSO (MEE  )  )  MEE-MEE* 1 

READ  (  I2)NN3,(  (S(LI  ,L2)  ,L1»1  ,IDIM) ,L2-I  ,NN3) 

Read  (  1  n  )  N  N 1  ,(  (S(LI  ,  L  2  )  *  L  l  •  1  ,  I  D  I  M  )  ,  L  2«  I  ,NN1  )  , 

I  I  I  B  I  U  I  ,1.2)  ,Ll*l  •  I  D  I  M  )  ,  L  2  •  I  ,  I  D  I  M  )  ,DET 
READ  I  1 1  )NN1  , ( ( A(LJ ,L2)  i L  1  *  1  ,  IDIM) ,  L2  a  1  »NN|  ) 

NN-NNI/101M 
NSPT  (  I  )  .NROW (  I  ) ♦NN 
NRWaNROW (  I  } -  I 
00  5  J" I , NN 
NRWaNRW* | 

IF  ILOC(NRIN)  .LT,  MISTY1G0  TO  5 
NSPT  (  1  )  » N R W 

NNO»(NRw-NROw(  I  )  )  •  I  0  I  M  ♦  I 
DO  B  L 1 ■ 1  , 1 0  1  M 
00  B  L2»NN0  ,  NN ] 

B  S  (  L  I  ,  L  2  )  ■  A  (  L  1  ,  L  2  ) 

BACKSPACE  10 

WRITE  ( I  0 ) N  N  I  ,  (  (S(LI  iL2) , L 1 ■ I ,!OlM)  ,L2»1  »NNJ  )  , 

1  I  (B(LI  ,  L  2 )  ,  L  1  ■  I  , I D 1 M  )  , L  2  » I  ,IDIM),OET 
IF  (1  ,FQ,  MRX ) GO  TO  6 

backspace  10 

READ  ( 1 0 ) NN  I  ,  <  (  S(L 1 ,L2)  • L I  *  1  , IDIM) ,L2»1  ,NN1  )  , 

1I(B(L1,L2),LI"1,IDIM),L2"I,ID|M),DET 

GO  TO  6 

5  CONTINUE 

6  PRINT  1 0 1  ,  1  , NROW (  J  )  ,  NSPT  (  J  ) 

'°7  CONTINUe'  1  "  **  NR0*,(,)  ■’•IS.'  ••  NSPT(I)  -*,15) 

REWIND  10 

B  DO  39  I-MISTY.NTNN 
NSPT ( 1 ) -NROW ( I ) 

READ  <12)NN3,((S(LJ,L2),lI»1,IDIM),l2»I,NN3) 

NN»NN3/ I D I M- 1 

IF  (MEE  ,GT.  NSPD ) GO  TO  17 

IF  (NREO(I)  .NE.  NWSO ( MEE  )  ) GO  To  11 
DO  10  Kal , IDIM 

IF  (SPD(MEE.K)  ,GE.  IOO.)GO  TO  10 
DO  9  L ■ 1  , 1 0  I M 


IB) 


54. 

S(L,NN3)3S<L.NN3)-S<L.K)*SPD(MEE,k> 

57. 

S(K  ,L  1-0. 

58. 

9 

5  (  L  ,  K  1  ■  0  . 

59, 

S  (<,<)■  1  , 

60. 

S  (  X  ,  NN 3  I »5P0 ( M£E  ,K  ) 

. 

1  0 

CONTINUE 

62. 

MEE-MEE* 1 

63. 

IE  <ME£  .GT,  NSPD1G0  TO  17 

69. 

1  1 

NU-  1 

65. 

NRW 1 -NROW ( 1 1-1 

6  6, 

DO  16  k-mee.nspd 

67. 

NWS-NWSO  <K 1 

6  8. 

00  19  L-NU.NN 

69. 

NRW-NRW 1 *L 

70. 

ir  (LOC(NRW)  .GT.  NOER ( NWS ) ) GO  TO  16 

71  . 

IE  (  L  0  C  (  N  R  W )  . N  E  «  NOER ( NWS  )  1  GO  TO  19 

72. 

mm«l 

73. 

M  3  «L • I 0  I  M 

79. 

00  13  !  A*  1  ,  10  1M 

75. 

IE  (  5  P  D  (  K  ,  l A )  .GE,  100. 1G0  TO  13 

76. 

DO  12  in-| , 1 0 1 M 

77  . 

S < 18  ,NN3 »-S< IB.NN3 »-S< IB .M3*t A )*SPD(K . I  A ) 

78. 

IE  (  > ( 1 n , 1 8 ) -  1 »  «  G  T .  . 001  )S( I B  ,M3*I A)«0. 

79. 

I  2 

CONT I NUE 

80. 

1  3 

CONT 1 NUE 

8  1  . 

GO  TO  15 

82. 

I  9 

CONTINUE 

83. 

GO  TO  17 

89, 

1  5 

N  U  ■  M  M  ♦  1 

85. 

1  6 

CONT 1 NUE 

8  6, 

1  7 

IE  (I  ,EQ.  1 IGO  To  36 

87. 

JILL-  1 

88  . 

1  1  -1-1 

89. 

DO  39  Km  1  , i  1 

9n, 

READ  (10)NN1,((A(LI,L2).L1»I|101M).L2-1,NN1), 

91  . 

1  <  <B(L1  .  L  2  1  ,  L  l  ■  1  ilDIM)  ,  L  2  ■  l  ,  l  D I M ) . OE  T 

92, 

I  SPT-NSPT (K ) 

93. 

IE  (LOClISPTI  .NE.  I ) G  0  To  38 

99  . 

I  TER  -  ( NSPT  <K ) -NRO* ( K )  ) • I 0  I  M 

95. 

NSPT (K  )  ■ N  SP  T ( K ) ♦ I 

96. 

IE  (JILL  .GT.  NSPDJGO  TO  22 

97. 

CO  18  L  1  -JILL ,NSPD 

98. 

NWO-NWSO (L 1 ) 

99. 

KEE-L 1 

l  on. 

IE  (NOER(NWD)  .EO.  <)G0  To  21 

101  . 

IE  (  NOE  R ( NWD 1  . GT  ,  K)GO  TO  22 

102. 

1  8 

CONT I NUE 

103. 

2  1 

JILL-<EE*1 

109. 

00  20  I  A  -  1 . I D 1 M 

105, 

IE  <SPO(KEE.m  .GE.  100. )G0  TO  20 

106, 

00  19  I  H- 1 , I D I M 

107. 

S(  IH,NM3)-5( IB.NN3)-SP0(KEE,IA)-A(Ia,ITER*IB) 

108. 

1  9 

a  (  i  a  .  i  ter ♦  i  b ) -o • 

109, 

20 

CONTINUE 

1  in. 

C 

GAUSS  IAN  ELIMINATION 

111. 

22 

00  25  I  A ■ 1  . I D I M 

1  12. 

I 0- l TE R  + I  A 
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113. 

Ill, 

115. 

1  1  6  . 
117. 

ns. 

1 1’, 
l  2n . 
121  . 
122. 
123. 
129. 
125. 

1  26. 

I  27. 

1  2  R , 

1  29. 

1  30. 
131. 

1  32, 
133. 

1  3*4. 

1  35. 
13  6. 

1  37. 

1  3  R  . 
1  39. 

1  MO . 

m . 

1  92. 
193. 

1  9r, 

195. 

196. 

197. 
1  9fl, 
1  99, 
150. 
151  . 

152. 

153. 
159, 

155. 

156. 

157. 

15  6. 

159. 

160. 

16  1. 
1  62. 
163. 
169, 

165. 

166. 
167. 
1  6  A , 
169, 


00  25  I B« I  , I  0 1 M 
C  <  1  A ,  1 B ) *0 , 

00  23  1C-1 ,101M 

23  C<!A,lBl-C<lA,IBl*AUC,|D)*BnC,!B) 

25  Continue 
NNOsNN  1  -  I  TER 

00  27  1  A«1  ,  I  0  I  M 

DO  27  IB»1,NN0 
0  (  1  A  i  ! B ) »0  . 

1  o  ■  1  T  F,  R  *  1  B 
00  26  1  C  ■  1  ,  I  0 1 M 

26  0(  1  A , 1B)»D( 1  A  *  1 B 1 ♦ C ( t  A • 1 C ) • A  < !C, !0) 

0 (  ! A  ,  1 B ) »D ( 1  A , 1 B 1  /DET 

27  CONTINUE 

L  Y  N  ■  N  N  0  /  1  0  1  M  -  1 
IF  (LYN  ,EQ.  0 1  GO  TO  32 
NRW-NROW ( 1 1-1 
JS  ■  1 

00  31  L«1,LTN 
MSPT« 1 5PT>l 
L  3»L • 1 0 1 M 
DO  29  M.JS.NN 
NRW.MRW+l 

IP  ( U  OC ( MSP  T 1  .NE.  LOC(NRM)  1  GO  TO  29 

M3» ( NRW.NRO*  < 1  1 ♦ l ) • 10  I M 
00  2fi  1  A  -  1  ,  : D  1  M 
00  28  1  B ■ 1  ,  I  0  I M 

1  L  ■  L  3  ♦  1  R 
1 M»M3  *  1 B 

2  R  S  (  1  A  ,  1 M  1  «S <  I  A , 1 M 1 -0  (  1  A , !L 1 
GO  TO  30 

29  CONTINUE 
GO  TO  32 

30  JS-mm+1 

31  CONTINUE 

32  00  3  3  1  A  » 1  , 1 0 1 H 

S{  1A,NN3)»S<IA|NN3)-0(IA,NN0) 
no  33  1 B ■ 1 , 1 0 ! N 
1  D  *  1  T  E  R  +  1  B 

A  (  |  A  ,  1 0 ) »C  <  I B . 1  A  1 /OET 

33  S  (  1  A , IB  » - 5  < I  A , 1B1-D ( I  A,  IB) 

BACKSPACE  10 

WRITE  <10)NN1,((A(L1iL2).L1«1,IDIM»,L2"1.NN1), 
1 <  (B(L1  »L2)  ,  L 1 ■ 1 • I D 1 M |  ,L2"1  ,  1  0  !  H  )  ,OET 
IF  ( K  ,EQ.  I  1  1  GO  TO  39 
BACKSPACE  10 

READ  C lOJNNl . ( < A (LI .L21  »Ll*l  •  I  0 1 H  > ,L2«1 ,NN1 ) , 
1  I  <P<ll*l2),Ll"l»!D!M)  iL2“1,IOIM),oET 
39  CONTINUE 

36  DET»5<  1  ,  1  )  *S ( 2 , 2 1 -S ( 1  , 2 ) *S ( 1  *2) 

IF  <  I  0  1  M  , EQ •  3 1  GO  TO  69 
H  <  I  , 1 ) -S  <  2 ,2 ) 

B(2,2)-S(  1,1> 

B (  1  , 2 1 »-S  < 1,2) 

GO  TO  65 

69  H(1,1)-S<2,2)*S(3,3»-S(2.3)*S(2,3) 


..  .  •  — 
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l  7n. 
171  . 

172, 

173, 
17% 
1  75. 
1  76, 
177. 
I7fl, 
1  79. 
1  AQ . 
I  8|  , 
182. 
1  83. 
1  89. 
1  85, 
186. 
1  87  . 
188. 
1  89. 
1  90. 
1  9  1  » 
1  92. 
1  93. 
1  9% 
195. 
1  96. 
197. 
1  9fl, 

1  99, 
20n, 
201  , 
202. 
203. 
201 . 

205. 

206, 
207, 
20fl. 
209. 

2  10. 
2  11, 
212. 
213. 
219, 
2  15. 
216. 


8(  1  ,  2  )  «s ( l  ,  3  > • S ( 2 , 3 ) «S (  1  • 2  1  *S ( 3  ,  3 ) 

8(  1  .3>»S(1,2)*SC2,3)-S(1  |3»*S<2,2> 

8 ( 2 , 2  )  -  S ( 1  .1  ) *S ( 3 , 3 ) »S (  1 »  3 1 *S ( 1  ,3) 

B  (  2 .3 ) >s ( 1 , 2 1  * S (  1  . 3 ) *S (  1  i  1 ) *S ( 2 ,3 ) 

B ( 3 . 3  )  -0ET 

DETa5(l,3)*8(l,3)+S(2,3)*B(2»3)*S(3,3)*B(3.3) 

B ( 3 ,2  1  -R (2,3  > 

B  (  3  ,  1  )  «  B  (  1  ,  3  > 

65  B(2, 1  )«R( i  ,2) 

1 F  (I  .EQ.  NTNN ) 60  TO  39 
NN 1  a  N  N  •  ( 0 1 M* 1 
N  3  ■  1  D  I  M  ♦  1 

WRITE  (ll)NNl,((S(Ll,L2),LlPl,IDlM),L2aN3,NN3) 
WRITE  (l0)NNl,((S(Ll,L2),Ll*l»10lMl,L2aN3,NN3), 
1  (  (B<L1  ,1.  2)  ,  L  1  *  1  •  1  D  1  M  )  ,L2"1  , 1 D 1 M ) . DE  T 
rewind  in 

39  CONTINUE 

c  back-substitution 

0  0  76  1  a  1  , l D  1  M 
0 ( I  ,  1  )  -0. 

DO  7  1  J« 1  ,  I  D IM 

71  D( 1 , 1 1  -  D  ( 1 ,1 )  +B  ( 1 ,J)*S(J,NN3) 

76  0(  1  ,11-0  (1,1)  /OET 
DO  75  1-2.NTNN 
NODE  *NTNN- 1 ♦ 1 
BACKSPACE  10 

READ  ( 1 0 ) N  N 1  ,  ( (A(L1  » L  2  1 • L  1  ■  1  ,  1  D  1  M  ) ,  L  2  ■  1  ,NNl  )  , 

1  ((B(L1,|.  2),Ll-l  ,lDIM),i.2*l  ,1D1M),0ET 
backspace  io 

on  77  i a ■ I , 1 0 1 h 
0  (  I  A  ,  1  1  a  0  « 

00  72  1 B- 1  ,  1  D  1  M 

7?  0(  1  A  ,  1  )  *0 (  1  A ,  I  1  ♦  B  (  1  A |  I R  )  •  A (  1  B  ,  NN 1  ) 

77  D (  I  A  ,  1  1  -0( 1  A , 1 ) /DET 
N  R  W 1 -NROW ( NODE ) -  1 
NN«NN 1 / I o  1  M 

DO  79  J« 1  ,NN 
JJ- ( J- 1 ) • 1 0 1 M 
NRW-NRW 1  * J 
MRW.NTNN-LOC ( NRW ) ♦  1 
DO  73  1  A ■ 1  , I D I M 
DO  73  1B-1 , ID1M 

73  D{  I  A  .  1  )  -  D  (  1  A , 1  )-A(  1A.JJ«1b)*D( IB.MRW) 

79  CONTINUE 
75  CONTINUE 
return 
END 


end  or  COME  1 L AT  1  ON: 


NO  DIAGNOSTICS, 
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APPENDIX  B 

DATA  CARDS  PREPARATION  GUIDE 

The  sequential  arrangement  of  the  data  cards  supplied  during  each 
program  run  as  well  as  the  input  items  punched  on  each  card  are  shown  in 
Table  A.  Regarding  the  input  parameters,  the  following  explanatory  notes 
are  made: 

Card  II  1 

IS0  -  1  in  anisotropic  cases;  0  in  nonhomogeneous 

cases;  -1  in  isotropic  cases. 

NR  -  Total  number  of  annular  surfaces  dividing  the 

Brazilian  test  cylinder. 

NS  -  Total  number  of  sectors  the  circular  face  of 

the  cylinder  or  its  quadrant  is  divided  into. 

NZ  -  Total  number  of  divisions  in  the  Z-direction; 

0  in  2-aimensional  problems. 

NW0P  -  1  if  whole  circular  face  of  cylinder  is  involved 

in  analysis;  0  if  only  a  quadrant  is  involved  in 
analysis. 

NDIM  -  2  in  2-dimensional  cases;  3  in  3-dimensional 

cases . 

NCYC  -  Total  number  of  load  cycles. 

NER  -  Least  total  number  of  similar  and  consecutively 

numbered  elements  in  finite  element  mrsh;  1  in 
nonhomogeneous  and  anisotropic  cases.  This  para¬ 
meter  is  irserted  to  avoid  having  to  completely 
treat  each  element  all  over  again  even  if  all  are  the 
same . 

NSS  -  Difference  between  the  numbers  of  a  loaded  top 

node  and  the  diametrically  opposite  bottom  node; 

0  if  NW0P  is  0  or  if  finite  element  mesh  is  sym¬ 
metric  about  plane  of  diametral  loads. 


Table  A  Input  Cards  Arrangement 


Spaces  for  EMI,  EM2,  CS1,  CS2,  PR1,  PR2  are  left  blank  in  anisotropic  cases. 
Supplied  in  3-dimensional  cases  only. 

Supplied  in  nonhomogeneous  cases  only. 

Supplied  in  anisotropic  cases  only. 
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Card  II  2 

TCRT 

Ratio  of  allowable  tension  to  allowable  com¬ 
pression. 

SCRT 

Ratio  of  allowable  shear  to  allowable  com¬ 
pression. 

EMI,  FM2 

Range  of  values  of  elastic  moduli. 

CS1,  CS2 

Range  of  values  of  allowable  compressive 
stresses. 

PR1,  PR2 

Range  of  values  of  Poisson's  ratios. 

Card  H  3 

P0SM 

Factor  to  be  multiplied  to  the  stiffness  matrices 
of  failed  elements  to  obtain  that  portion  of  the 
stiffness  matrices  to  be  subtracted  from  the 
global  stiffness  matrix. 

TRN 

Maximum  deviation  which  an  element  load  factor 
can  have  from  the  critical  load  factor  for  the  ele¬ 
ment  to  be  considered  failed. 

DIA 

Diameter  of  cylinder. 

H 

Length  or  thickness  of  cylinder. 

NTEP 

0  if  failure  criterion  in  which  the  elastic  modulus 
across  tension  cracks  is  reduced  to  zero  is  to  be 
applied;  1  if  the  old  failure  criterion  (described 
in  first  annual  report)  is  to  be  applied;  always  1 
in  anisotropic  cases. 

NG0T 

0  if  both  finite  element  mesh  and  material  property 
are  symmetric  about  a  nodal  plane  perpendicular 
to  the  Z-axis;  1  if  no  symmetry  exists.  With  this 

parameter,  only  one-half  of  the  cylinder  lengtn  or 
thickness  need  be  analyzed  if  symmetry  exists . 
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NST1 ,  NST2 
MST1,  MST2 


Pairs  of  points  between  which  magnitudes  of 
strain,  usually  the  horizontal  strain  at  the 
center  point  of  each  end  of  the  cylinder,  are 
desired . 


Card  H  5a  ,  etc . 

RC(1)  -  0  for  solid  cylinders;  inner  radius  of  hollow 

cylinders . 

RC(2)-»  RC(NRl)  -  Radial  coordinates  of  annular  surfaces. 


Card  #  6a,  etc. 

TC(I),  1=1,  NS1  -  Circumferential  coordinates,  in  degrees,  of 

nodal  radial  planes  measured  clockwise  from 
plane  of  loads.  TC(1)  =  0. 


Card  II  7a.  etc . 


ZC(I),  1=1,  NZ1 


Z-coordinates  of  nodal  circles  using  an  end 
circle  as  datum.  ZC ( 1)  =  0. 


Card  H  8 

N  -  Any  10-digit  integer.  This  is  the  starting  point  of 

random  number  generator. 


ALF2 ,  ALF3 


Angle,  in  degrees,  which  the  line  A'B'  (see 
Figure  2.4  of  semi-annual  report)  makes  with  the 
plane  of  the  loads  . 

The  angles  and  0:3,  in  degrees,  defined  in 
Figure  2.4  of  semi-annual  report. 

These  three  angles  define  the  directions  of 
anisotropy. 
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Cl,  C2 ,  C3  -  Allowable  compressive  stresses  in  the  direction 

of  axes  of  anisotropy. 

Cl,  E 2  -  Elastic  moduli  in  the  direction  of  axes  of 

anisotropy,  1,2. 


Card  It  10 


C3 


PI,  P2,  P3 
G 1 ,  G2 ,  03 


Elastic  modulus  in  the  direction  of  axis  of 
anistropy  3. 

Poisson's  ratios  along  axes  of  anisotropy. 

Shear  moduli  associated  with  directions  of 
anisotropy. 


The  entries  below  the  heading  F0RMAT  in  Table  A  are  called  format 
specifications.  These  serve  as  instructions  to  the  user  on  how  the  input  data 
are  to  be  punched  in  the  data  cards.  In  the  present  program,  each  item  in  the 
input  list  has  the  format  specification  of  the  form  Fw.d  or  Iw  in  which  w  is 
the  number  of  column  spaces  on  the  card  reserved  for  each  item.  Given  the 
w  s  of  the  entire  input  list,  it  is  an  easy  matter  to  determine  the  exact  columns 
on  the  card  where  each  input  item  is  to  be  punched.  Tor  example,  in  card  ti  3 
one  punches  P0SM  in  columns  1-10,  TRN  in  columns  11-20,  DIA  in  columns 
21-30,  H  in  columns  31-40,  NTEP  in  columns  4  1-45,  and  NG0T  in  columns 
46-50.  Integers  (items  with  I  specifications)  are  punched  right-justified  with¬ 
in  the  spaces  allocated  to  them.  Floating-point  numbers  (items  with  F  specifi¬ 
cations)  may  be  punched  with  or  without  the  decimal  point  anywhere  within 
their  allocated  spaces.  If  punched  without  a  decimal  point,  the  computer  will 
read  the  floating-point  number  as  if  there  was  a  decimal  point  d  places  to  the 
left  of  the  right-most  digit.  The  number  before  a  F  or  an  I  is  a  repetition  factor; 
for  instance,  (415)  =  (15,15,15,15). 

Input  data  listings  are  shown  in  Chapter  2. 


